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Abstract 

Based on an approximate six-quark operator effective Hamiltonian from perturbative QCD, 
we present a systematical study of charmless B PP, PV, VV decays {P and V denot- 
ing pseudoscalar and vector mesons, respectively). The calculation of the relevant hard- 
scattering kernels is completed, the resulting transition form factors are consistent with 
the results of QCD sum rule calculation. Important classes of power corrections include 
"chirally-enhanced" terms, vertex corrections and weak annihilation contributions with non- 
trivial strong phase. With these considerations, predictions are presented for the branching 
ratios and CP asymmetries of B-meson decays into PP, PV and VV final states, and also 
for the corresponding polarization observables in VV final states. Several decay modes and 
observables, which are of particular interest phenomenologically, are discussed in detail, in- 
cluding the effect of annihilation amplitude with strong phase, the tttt, ttK and irp systems, 
the longitudinal polarization fraction in pK* and (j)K* systems and so on. It is observed 
that predictions in our framework generally agree with the current experimental data. 



1 Introduction 



The study of hadronic charmless B-meson decays can provide not only an interesting avenue to 
understand the CP violation and the flavor mixing of quark sector in the standard model (SM), 
but also a powerful means to probe different new physics scenarios beyond the SM [HE]. With the 
operation of dedicated B-factory experiments, a huge amount of experimental data on hadronic 
B-meson decays has been analyzed with appreciative precision. To account for the experimental 
data, theorists are urged to gain deeper insight into the mechanism of these decays. 

Theoretically, in order to consistently predict hadronic B decays, it needs to deal with 
the short-distance contributions in a complete and systematic way from the high energy scale 
down to a proper low energy scale at which the perturbative calculations still remain reliable, 
and to treat the long-distance contributions which contain the non-perturbative strong inter- 
actions involved in these decays. The main task is to reliably compute the hadronic matrix 
elements between the initial and final hadronic states. In the past years, much progress has 
been made in our understandings of the hadronic charmless B-meson decays Pj, and several 
novel methods based on the "naive" factorization approach (NF) [3j, such as the perturbative 
QCD method (PQCD) [1], the QCD factorization approach (QCDF) and the soft-collinear 
effective theory (SCET) [B], have been proposed. 

All the above frameworks of weak decays are based on the four-fermion operator effective 
Hamiltonian derived via operator product expansion and renormalization group evolution. In 
hadronic weak decays, the short-distance QCD contributions are characterized by the Wilson 
coefficient functions of four-quark operators and the long-distance ones are in principle obtained 
by evaluating the hadronic matrix elements of relevant four-quark operators. The Wilson coeffi- 
cient functions are in general calculated by perturbative QCD which is well developed, while the 
evaluation of hadronic matrix elements remains a hard task as it involves non-perturbative effects 
of QCD. In fact, for the mesonic two-body decays of B meson, it involves three quark-antiquark 
pairs once each meson is regarded as a quark-antiquark bound state at the quark-level structure. 
This fact then naturally motivates us to consider the six-quark (rather than four-quark) operator 
effective Hamiltonian [7J. 

For the infrared singularity caused by the gluon exchanging interaction during the evaluation 
of the hadronic matrix elements of effective six-quark operators, it is simply regulated by the 
introduction of a mass scale motivated from the gauge invariant loop regularization method [8] , 
where the energy scale fig is introduced to play the role of infrared cut-off energy scale without 
violating gauge invar iance. 

We have applied the QCD factorization based on six-quark operator effective Hamiltonian 
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to —7- tttTjItK, KK decays [7J, and the theoretical predictions for all the branching ratios 
and CP asymmetries in these decays are found to be consistent with the current experimental 
data except for a few decay modes. The strength of annihilation diagram contributions is a 
widely discussed issue in hadronic B decays [9]. Their effects are assumed to be important both 
in pQCD and QCDF framework, and could even be fine-tuned independently for different decay 
modes \10\ lllj . Motivated by these observations, in this work we try to find a simultaneous 
solution for all B — )• PP, PV, VV decay modes in the framework of six-quark operator effective 
Hamiltonian, by considering annihilation contribution with non-trivial strong phase. However, 
the calculation of strong phase from nonperturbative QCD effects is a hard task, there exist no 
efficient approaches to evaluate reliably the strong phases caused from nonperturbative QCD 
effects , and we has to set the strong phase as an input parameter in our framework. 

Our paper is organized as follows. In section[2l after briefly reviewing the four-quark operator 
effective Hamiltonian, we begin with the introduction of the primary six-quark diagrams with 
the exchanges of a single W-boson and a single gluon, and the corresponding initial six-quark 
operator. It is shown that a complete six-quark operator effective Hamiltonian is in general 
necessary to include all contributions from both perturbative and non-perturbative QCD cor- 
rections. The treatments of the singularities caused by the gluon exchanging interactions and 
the on mass-shell fermion propagator are presented in Section [3l Then in the next section, the 
vertex corrections and annihilation contributions are presented. Section [5] contains all the input 
parameters which will be used in our calculation. In Sectional we will give our numerical predic- 
tions and discussions for B — )• PP, PV, VV decays. Our conclusions and the decay amplitudes 
are given in the last section and in the Appendix, respectively. 

2 Effective Hamiltonian of Six-Quark Operators 
2.1 Four-Quark Operator Effective Hamiltonian 

Let us start from the four-quark effective operators in the effective weak Hamiltonian. The 
initial four-quark operator due to weak interaction via W-boson exchange is given as follows for 
B-meson decays 

Oi = {qibi)v-A{^jUj)v-A, q"" = u, c, q"^ = d, s. (1) 

The complete set of four-quark operators are obtained from QCD and QED corrections 
which contain the gluon-exchange diagrams, strong penguin diagrams and electroweak penguin 
diagrams. The resulting effective Hamiltonian (for 6 — )• s transition) with four-quark operators 
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is given as follows [12j 

r \ 10 ■ 



* q=u,c 
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where = Vgh^qs a-re products of the CKM matrix elements, Ci{^) the Wilson coefficient 
functions jl2]. and Oi{^) the four-quark operators 

= {sibi)v~A{qjqj)v-A , 

'^{qA)v-A(.Sjqj)v-A , 
q' 

-- -'^y2i^'A)s-pisjqj)s+p , 

V (3) 

= -^^eg,{qibi)s-p{sjq'j)s+P , 

q' 

q' 

Here the Fermi constant Gp = 1.16639 x 10~^ GeV~^, {q'q')v±A = 7^1(1 iTs)*?') and i,j are 
the color indices. The index q' in the summation of the above operators runs through u, d, s, 
c, and b. The effective Hamiltonian for the b ^ d transition can be obtained by changing s into 
d in Eqs. ©and ([3]). 

2.2 Six-quark Diagrams and Effective Operators 

As mesons are regarded as quark and anti-quark bound states, the mesonic two body decays ac- 
tually involve three quark-antiquark pairs. It is then natural to consider the six-quark Feynman 
diagrams which lead to three effective quark-antiquark currents. The initial six-quark diagrams 
of weak decays contain one W-boson exchange and one gluon exchange, thus there are four 
different diagrams as the gluon exchange interaction can occur for each of four quarks in the 
W-boson exchange diagram, see Fig. [TJ 

. h h 



Figure 1: Four different six-quark diagrams with a single W-boson exchange and a single gluon 
exchange. 

The resulting initial effective operators contain four terms corresponding to the four dia- 
grams, respectively. In a good approximation, the four quarks via W-boson exchange can be 
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regarded as a local four-quark interaction at the energy scale much below the W-boson mass, 
while the two QCD vertexes due to gluon exchange are at two independent space-time points, the 
resulting effective six-quark operators are hence in general non-local. The six-quark operators 
corresponding to the four diagrams in Fig. [3] are found to be 

0(f = 4vra, ///^-^^e-«---^)P+(---3W(^-,(^3)^^rV(^^ ^ 



— m? + ie 



xfe(:r2) 7"r"rigi(xi)) * fe(xi)r2(?3(^i)), 

p'^ — rrig^ + i€ 

0(6) = 4™, //^-^e-^«-•l--■2)P+(^2-^3)'=)(g'(x3)7.TV(a:3))- 



— -|- le 

xfe(rci)ri,i(li)) . fefe) Z'^"'"' . 7"rT293(rCi)). (4) 

— m^^ -|- le 

where k and p correspond to the momenta of gluon and quark in their propagators, qi is usually 
set to be the heavy b quark, xi, X2 and X3 are space-time points corresponding to the three 
vertexes. The color index is summed between ^1,52 and 93,(74. Note that all the six-quark 
operators are proportional to the QCD coupling constant Us due to gluon exchange. Thus the 
initial six-quark operator is given by summing over the above four operators 

o(6)=x:og). (5) 

Actually, the initial six-quark operators {j = 1,2,3,4) can be obtained from the following 
initial four-quark operator via a single gluon exchange 

o = ferigi) * (g4r2Q'3)- (6) 

Unlike the classical four-quark effective operator, the six-quark operators used here are non- 
local with quark and gluon propagators inserted into them. Such operator is equal to the one 
used in SCET with and be found to be a more effective form to describe dynamics at low 
energy scale around 1.5 GeV. 

With the above considerations, the QCD factorization approach with six-quark operator 
effective Hamiltonian enables us to evaluate all the hadronic matrix elements of nonleptonic two- 
body B-meson decays. The detailed calculations of the hadronic matrix elements in i? — )• 
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as an example, could be found in our previous paper [7j. As for the hadronic matrix elements 
in other B — )• PP, PV, VV decays, we shall list them in the Appendix. 



3 Treatment of Singularities 

Before proceeding, we would like to point out that there are two kinds of singularities in the 
evaluation of hadronic matrix elements. One singularity stems from the infrared divergence 
of gluon exchanging interaction, and the other one arises from the on mass-shell divergence of 
internal quark propagator. 

In general, a Feynman diagram will yield an imaginary part for the decay amplitudes when 
the virtual particles in the diagram become on mass-shell, and the resulting diagram can be 
considered as a genuine physical process. It is well-known that when applying the Cutkosky 
rule il3j to deal with a physical-region singularity of all propagators, the following formula 
holds: 

^ \ I . = P 2 -^^'^[Z - "T-q]; (7) 

— -|- Yp'^ — ^ 

which is known as the principal integration method, and the integration with the notation of 
capital letter P is the so-called principal integration. 

However, the Cutkosky rule is useless for singularities from infrared divergence of gluon 
propagator. Integration with those propagators is sensitive to the infrared cut-off for gluon and 
light-quark propagator, and diverge to infinity when the cut-off becomes to zero. A modified 
integration with different parameters for different channels is used in QCDF framework [lOllllj. 
while the transverse momentum kj- dominating in the zero momentum fraction is added to the 
propagator in pQCD framework [3]. In this work, we prefer to add the same dynamics mass for 
both gluon and light quark to investigate the infrared cut-off dependence of perturbative theory 
prediction: 

"PJT"^ ^ 77^5 o 5 ^ (q IS a light quark). (8) 

It is noted that, as the gauge dependent term kf^kiy can always be written as linear combina- 
tions of the momenta pa on the external lines of spectator quark, and they are all on mass-shell 
in our present consideration (as defined in Fig. [21 their contributions equal to zero once the 
equation of motion is used. Our results are therefore gauge independent. 
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Figure 2: Definition of momentum in — >• M1M2 decay. The light-cone coordinate is adopted 
with (n+, n^, k±). 

4 Vertex Corrections and Annihilation Contributions 

The short-distance contributions characterized by the Wilson coefficient functions for the effec- 
tive four-quark operators were calculated by several groups at the leading order (LO) and the 
next-to-leading order (NLO) [12l[Tl]. Their values depend mainly on the choice for the running 
scale /i. When considering the NLO Wilson coefficient functions and a^, one needs to include 
the magnetic penguin-like operator Osg which is defined as [12] 

Osg = ^mi,q,a^,{l + ^^)T^^G''^'^hj , (9) 

where i, j are the color indices. The magnetic-penguin contribution to a generic B — t- M1M2 
decay leads to the modification for the Wilson coefficients corresponding to the QCD penguin 
operators, 

^ «4,6(/i)-^^Qt(/.), (10) 



with = Csg + C5, \P\ = m'^/2, and 04^5 = ^^4,6 + which appear in the factorizable 
diagrams. 

As shown in Ref. CP violating observables may be improved by adding vertex correc- 
tions. Furthermore, the vertex corrections were proposed to improve the scale dependence of 
Wilson coefficient functions of factorizable emission amplitudes in QCDF [16] . Those coefficients 

ana L^2n + tvcT 

corrections, are modified to 



are always combined as 6*2^-1 + jj^ and C2n + '^'m ^ , which, after taken into account the vertex 



^^2n-l(/u) + ^(/Z) ^ C2„_l(/Z) + ^(/Z) + ^CF^^V2n-,{M2) , 

CM + ^(M) ^ C2M + ^(M) + ^C^^^V2.(M2) , (11) 

with n = 1, 5, M2 being the meson emitted from the weak vertex. In the naive dimensional 
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regulation (NDR) scheme, Vi{M) are given by [TOl [T6] 
Vi{M) = 



' 121n(^)-18 + /J(ix,/)a(x)5(2;) , for i = 1 - 4, 9, 10 , 

_121n(2^)+6-/o , fori = 5, 7, (12) 

—6 + Jq^ dx <pb{x) h{x) , for i = 6, 8 , 



where (pai^) and </>;,(x) denote the leading- twist and twist-3 distribution amplitudes for a pseu- 
doscalar meson or a longitudinally polarized vector meson, respectively. While for a perpendic- 
ularly polarized vector final state, (pa{x) = 4'±{x,fi) and </>f,(x) = 0. The functions g{x) and h{x) 
used in the integration are given respectively as [LQl 



gix) = 3( ^Inx — ivr 

\ — x 



2 Li2(x) - \r? x + - (3 2i vr) In X - (x O 1 

1 — X 



+ 

,2 



= 2Li2(x) -In^x - (1 -F2i7r)lnx- (x o 1 -x) . (13) 

To further improve our predictions, we shall examine an interesting case that vertexes receive 
additional large non-perturbative contributions, namely the Wilson coefficients ai = Ci + 
are modified to be the following effective ones: 

a, ^ = Q(^) + ^(/^) + ^Cf^^^{V,{M,) + V,{M,)) , (i = 1 - 4, 9, 10), 

a, ^ al^f = Q(^) + ^(^) + ^C^%J^(^.(M2) + y2(M2)) , (z = 5 - 8). (14) 

The corrections Vi(M2) and V2(M2) depend on whether the meson M2 is a pseudoscalar or a 
vector. It could be caused from the higher order non-perturbative non-local effects, as shown 
in FigEl Adopting Vi{P) = 26e-¥, V2iP) = -26, ViiV) = I5e¥, and V2iV) = -15et^ both 
the branching ratios and CP asymmetries of most B — t- PP, PV, VV decay modes are improved, 
which will be detailed later. 

After involving the vertex corrections and the contributions of effective Wilson coefficients, 
we next take account of the annihilation contributions with non-trivial strong phase. Both 
the vertex corrections and effective Wilson coefficients only affect emission diagrams while pall 
on exchange and annihilation diagrams. In our calculation, vertex correction for annihilation 
diagram might be sizable in PV and/or VV final state, and an extra strong phase for PV and/or 
VV channel would improve our theoretical predictions. 

Most of the annihilation contributions are from factorizable annihilation diagrams with the 
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Figure 3: The diagrams in (a), (b) and (c) are loop contributions only to the effective weak vertex (type 
I), only to the gluon vertex (type II), and for both weak and strong vertexes (type III), respectively. 



{S — P) X {S + P) effective four-quark vertex: 



4^^^(M). 
4p'(M). 



dxdy 
dxdy 
dxdy 
dxdy 



ifj-Pi + fip^)y{l - y) 


(x(l- 




- Hl + ie){{l -y)m% - ml 

- 3(2x - l)mv2)2/(l - y) 


+ tey 


(x(l- 


y)m\ 
(-3(1 


- ^il + ie){{l -y)m% - ml 

- 2x)mvi - /iP2)y(i - y) 


+ iey 


(x(l- 


y)m\ 


- fil + ie){{l -y)m\ - ml 


+ iey 


, 3(1 


-2x)(. 


-mvi + 3(23; — 1)771^2)7/(1 - 


-y) 


(xil- 


- y)m\ 


- ^2 + ie)((i _ y)m\ - ml 





(15) 



From Eq. (|15p and Eqs. (jA-4ip -( |A-43p . it can be ^ can be converted into 



A^p^a^^^ by exchanging x -H- (1 — y). Since the contributions of these amplitudes are dominated 
by the area x ~ or y ~ 1, A^^^-'{M) and Af},^'(M) have the same sign, while A^^f^{M) and 
A'i^p'^{M) have a different sign from A^'p^{M). As a result, we use the same strong phase for 
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A^]^%M) and A^'^^M), and another one for Aff^iM) and Af^%M). 

5 Theoretical Input Parameters 

The theoretical predictions in our calculations depend on many input parameters, such as the 
Wilson coefficients, the CKM matrix elements, the hadronic parameters, and so on. Here we 
present all the relevant input parameters as follows. 

5.1 Light-Cone Distribution Amplitudes 

For the wave function of B meson, we take the following form in our numerical calculations [17j: 



hix) = Nbx'^{1 - xfexp 



1 / xm^ Y 
'2 \ W) 



(16) 



where the shape parameter of B meson is tJB=0.30GeV, and Nb is a normalization constant. 

We next specify the light-cone distribution amplitudes (LCDAs) for pseudoscalar and vector 
mesons, the definition of which could be found in Refs. I19| . The general expressions of 
twist-2 LCDAs are 

oo 

Mx,l^) = 6x(l-x)[l + 5^a^(;u)C3/2(2x-l)], 

n=l 

oo 

^v{x,f^) = 6x{l-x)[l + ^a^{i^)Cl/\2x-l)], 

n=l 

oo 

</.^(x,//) = 6x(l-x)[l + ^a^'^(/i)C3/2(2x-l)], (17) 



n=l 



and those of twist-3 ones are 



0p(x,^) = 1, (p„{x,iJ,) = Gx{l - x), 

oo 

Mx,^^) = 3[2x-l + ^a^'^(/i)P„+i(2x-l)] 

n=l 

(f)+{x) = 3{l-xf, (/>_(x) = 3x^ (18) 



where C„(x) and -Pri(x) are the Gegenbauer and Legendre polynomials, respectively. The shape 
parameters of light mesons are taken from [2U] and listed in Table [TJ 

The parameters in Table[T]are given at the scale /x=1.0GeV, and should run to our calculation 
scale n = y/2AQCDmi, ~ 1.5 ± O.lGeV, with Aqcd ^ 288l^8MeV. It is noted that LCDAs of 
light mesons become more close to their asymptotic forms (all shape parameters become zero) 
when the scale runs to higher values. 
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Table 1: Values of Gegenbauer moments at the scale /i=lGeV, which are taken from |20) . 





TT 


K 


P 


K* 


<^ 




ai 




0.06 ± 0.03 




0.03 ±0.02 








0.25 ±0.15 


0.25 ±0.15 


0.15 ±0.07 


0.11 ±0.09 


0.15 ±0.07 


0.18 ±0.08 


T 








0.04 ±0.03 






T 
0^2 






0.14 ±0.06 


0.10 ±0.08 


0.16 ±0.06 


0.14 ±0.07 



5.2 Decay Constants and Other Input Parameters 

For decay constants of various mesons and other hadronic parameters, we list them in Table [21 
As for the CKM matrix elements, we shall use the Wolfenstein parametrization [25] with the 
four parameters chosen as [26l: A = 0.798t[];[]^^, A = 0.2252to:ooo82, p = 0.14ll°;°^f, and 
fj = 0.340 ± 0.016. 

Table 2: The hadronic input parameters [21] and the decay constants taken from the QCD sum 
rules [221 [23] and Lattice theory [21]. 







TUB 


mi, 






rrid 


1.638ps 


1.525ps 


5.28GeV 


4.4GeV 


173.3GeV 


4.2MeV 


7.6MeV 


rric 


rus 


m^± 




rriK 


nipO 


mp± 


1.5GeV 


0.122GeV 


O.UOGeV 


0.135GeV 


0.494GeV 


0.775GeV 


0.775GeV 










P-n 


PK 


f<f> 


1.7GeV 


1.8GeV 


300MeV 


0.78GeV 


1.02GeV 


0.892GeV 


0.215GeV 


fB 


U 


fx 


fp 




fK* 




0.210GeV 


O.lSOGeV 


0.16GeV 


0.216GeV 


0.187GeV 


0.220GeV 


0.151GeV 


Ik* 


J 4, 


Jp 










0.185GeV 


0.186GeV 


0.165GeV 











In our numerical calculations, the running scale is taken to be 

^l = ^2AQCDmb ^ 1.5 ± O.lGeV. 



(19) 



The scale of in the six-quark operator effective Hamiltonian is also taken at 1.5 GeV. Mass 
of b quark running to /x = 1.5 ±0.1 GeV follows the framework in |27] as: 



'Q^\7o//3o 

vr / 



TT 



27r2 



vr-' V 6 



(20) 



The definition of Cj can be found in [27J. After calculation, we can get mf,(/u) ~5.54 GeV. 
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In addition, the infrared cut-offs for gluon and light quarks are the basic scale to determine 
annihilation diagram contributions (the smaller /x^, the larger contributions), and are set to be 
fiq = /ig=0.37 GeV. 



5.3 Form Factors 

As is known, form factors have to be provided from outside in QCDF and SCET frame- 
works (such as by resorting to QCD sum rules or lattice QCD), while can be calculated in 
pQCD approach. The method developed here allows us to calculate the relevant transition form 
factors as below: 



with: 







V 



aB^Mi 
^0 



47ra(/i)C, 



F rpFNhhh, 



Ncm%FM2 



B){Mi,M2 = P), 



^M2 ("T-l 



Ncm%FM2 
47ra(/i)C 



mMimM2) 
V,M2 = P), 



-{Mi,M2 = V), 



F_rpFMiM2 



Ncmj^FM2 



LL,// 



m 



B 



-iMi,M2 = V), 



(21) 



'-LL,. 



Tll, 



Cf 



2Nr 



Before giving predictions of the observables in i? — )• PP, PV, VV decays, we first list our numer- 
ical results for the form factor at = 0, and for comparison, we also list the results of QCD 
sum-rules, light-cone sum rules and PQCD [281 [291 [30] in Table [3l 



6 Numerical Results and Discussions 

In this section, we shall classify the 43 channels of B decays into two light mesons according to 
the final states, and give our predictions for the branching ratios, the CP asymmetries, and the 
longitudinal polarization fractions. Comparisons with the current experiment data and other 
theoretical predictions, if possible, are also made. According to different decay modes, we shall 
give our predictions for the observables one by one. 

6.1 B ^ PP decays 

The contributions of effective Wilson coefficients and effects of different strong phases for anni- 
hilation diagrams for B — t- PP are shown in Tables [D and [5] For B — )• vrvr decay channel, the 
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Table 3: The B ^ P,V form factors at = in QCD Sum Rules, Light Cone and our work, where 
the errors stem mainly from the uncertainties in the global parameters fj.scaie — 1-5 ± O.lGeV, /Xg — 
0.37 ± 0.037GeV, and from the shape parameters of light mesons. Within their respective uncertainties, 
our predictions are consistent with the results of the other two methods. 



Mode 


F(0) 


QCDSR[30J 


LCSR[28J 


PQCD[29J 


This work 


B ^ K* 


V 


0.411 


0.339 


0.406 


n 977+0.075+0.008 
' -0.043-0.006 




0.374 


0.280 


0.455 


n 090+0.095+0.021 


Ai 


0.292 


0.274 


0.30 


n 990+0-054+0.006 

U.ZZU_Q Q3]^_Q QQ5 


B ^ p 


V 


0.323 


0.298 


0.318 


n 900+0.063+0.006 

U.ZOO_Q Q3g_Q QQ2 


Aq 


0.303 


0.248 


0.366 


n 9an+0. 077+0. 012 
'-'•^"'-'-0.041-0.005 


Ai 


0.242 


0.239 


0.25 


n 1 Qq+0.048+0.004 
U.iyo_Q Q28_o 002 


B ^uj 


V 


0.293 


0.275 


0.305 


n 9n«+0.056+0.005 
U.zuo_Q Q32_o.002 


Aq 


0.281 


0.231 


0.347 


n 9c:i +0.070+0.010 

U.ZJi_Q Q3g_Q 


Ai 


0.219 


0.221 


0.30 


n 1 70+0.043+0.004 
'-'•-'^''-'-0.024-0.002 


B -^ir 


Fo 


0.258 


0.285 


0.292 


n 9RQ+0.053+0.008 

U.ZUy_Q Q3Q_Q QQ5 


B^K 


Fo 


0.331 


0.345 


0.321 


n q/iQ+0.070+0.012 
'-'•'3^^-0.038-0.007 



naive power-counting based on factorization theory predicts: 

Br{7r-iT+) > Br{iT--K^) > Sr(7r°7r°). (22) 

However, from the current experimental data [31j, we can see that Br(7r^-7r+)~Br(7r^7r^) and 
Br(7r''7r'^) is much larger than theory prophecy. From Table [5l we can see that: with only 
the vertex and NLO corrections included, the branching ratio oi B ^ vr+vr^ and the ones of 
B^ — ?• vr+vr'^, vr'^vr'^ modes are bigger and smaller than the current experiment data, respectively; 
after taking into account the contributions of effective Wilson coefficients, which can reduce the 
former but enhance the latter two, our predictions are well consistent with the data. 

The decay B — )■ irK are dominated by penguin contributions. For the ratios defined by [32] 

^ 2Br{B+ ^ 7r^K+) ^ BrjB'^ ^ ^'K+j 

" Br{B+ ^TT+K^) ' " 2Br{B0 ^ttORO)- ^ ' 

There is a relation Rc = Rn ~ 1 if other diagram amplitudes such as annihilation are negligible 
compared with penguin emission amplitude. The current experimental data are Rc = 1.12ib0.07 
and Rn = 0.99 it 0.07. In our framework, we have Rc = 1.15 and i?„ = 1.13, which are in perfect 
agreement with the data. 

If we only consider the NLO contribution and vertex corrections, the direct CP violations 
for B^ — ?• ir^ K"^ , Tr~^ are wrong in signs compared to the experiments. To improve the 
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Table 4: The branching ratios (in units of 10~^) and direct CP asymmetries in i? ttK decays. 
The central values are obtained at /ig = /Xg=0.37GeV, the first error stems from the uncertainties in 
the global parameters ^JLscale = 1-5 ± O.lGeV, /Zg = 0.37 ± 0.037GeV, and the second from the shape 
parameters of light mesons. NLO'' and NL0''''''''(6'") stand for resuhs with "NLO correction+effective 
Wilson coefficients" and "NLO correction+effective Wilson coefficicnts+annihilation with strong phase" , 
respectively. 



Mode 


Data[3T1 


This work 






NLO+Vertex 


NLO'= 




NLO'^-''-''(5°) 


NLO'=-''-''(20°) 


B+ ^ TT+A'O 


23.1 ± 1.0 


22.5 


21.4 


19.0 


22.6 


25.9 


B+ ^ 7r"A'+ 


12.9 ±0.6 


12.8 


12.5 


11.2 


13.1 


14.9 


BO -^tt-K+ 


19.4 ± 0.6 


19.2 


19.5 


17.4 


20.5 


23.3 


B" ^ 7r''A0 


9.8 ± 0.6 


8.3 


8.4 


7.4 


8.9 


10.2 


Acp{n+KO) 


0.009 ± 0.025 


-0.006 


-0.006 


-0.006 


-0.007 


-0.007 




0.050 ± 0.025 


-0.053 


0.012 


0.003 


0.018 


0.034 


Acp{Tr-K+) 


-0.098 ±0.012 


-0.118 


-0.139 


-0.158 


-0.131 


-0.105 


Acp{7t0k0) 


-0.01 ± 0.10 


-0.052 


-0.139 


-0.143 


-0.138 


-0.137 




0.58 ±0.17 


0.699 


0.760 


0.768 


0.756 


0.745 



Table 5: The same as Table S] but for B Tnr,KK decay modes. 



Mode 


Data[3T] 


This work 






NLO+Vertex 


NLO<=^'^' 


NLO^-''-''(-40°) 


NLO=-^-^(5°) 


NLO'=^^(50°) 


B° -S> TT-1T + 


5.16 ±0.22 


7.1 


6.5 


6.00 


6.6 


7.6 


B+ ^ tt+ttO 


5.59 ±0.40 


4.1 


5.5 


5.51 


5.5 


5.5 




1.55 ±0.19 


0.3 


1.0 


1.11 


1.0 


1.0 


B+ -s> K+K° 


1.36 ±0.28 


1.7 


1.6 


1.0 


1.7 


2.2 


BO ^ rORO 


0.96 ±0.20 


1.5 


1.4 


0.7 


1.5 


2.2 


B° -s> a:+a— 


0.15 ±0.10 


0.09 


0.09 


0.09 


0.09 


0.09 




0.38 ± 0.06 


0.206 


0.266 


0.239 


0.260 


0.141 


Acp(7r+7r°) 


0.06 ±0.05 


-0.000 


-0.001 


-0.001 


-0.001 


-0.001 


AcpCttOttO) 


0.43 ± 0.25 


0.382 


0.453 


0.272 


0.485 


0.789 




-0.61 ±0.08 


-0.504 


-0.506 


-0.353 


-0.524 


-0.638 


Acp{K+KO) 


0.12 ±0.17 


0.101 


0.098 


0.041 


0.101 


0.106 


AcpiKORO) 


-0.58 ±0.7 


0.000 


0.000 


0.000 


0.000 


0.000 


Acp(K+K-) 




-0.184 


-0.184 


-0.184 


-0.184 


-0.184 
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predictions, we apply the effective Wilson coefficients to get more reasonable results, and add 
a strong phase = 5° for annihilation diagrams. From Table [H it is noted that effective 
Wilson coefficients and annihilation contributions with strong phase could bring the signs of 
Acp{t^^ K^) back to the right track, but do nothing with the sign of Acp{Tr~^K^)- 

As for tree diagram dominated channels B — )• vrvr, KK, the contributions of strong phase are 
ignorable since they are doubly suppressed by CKM and small value of this strong phase. If we 
want to consider effect of strong phase, a large one such as —40° or 50° are needed here. We list 
their effects in Table [SJ 

In order to better test our theoretical framework, we list the most recent predictions based on 
QCDF with strong phase effects [11] and the predictions from pQCD [15j approach in Tables [6] 
and [71 The first theoretical error in our calculations is referred to the global parameters of 
running energy scale n scale and the infrared energy scale /i^, and the second one is from the 
shape parameters of light mesons. 



Table 6: Comparisons of predictions between our framework and QCDF, pQCD methods in B tit:, ttK 
decays. 



Mode 


Data 31 


QCDF 11^ 


pQCD 15 


This work 








LO 


NLO+Vertcx 


LO 


NLO+Vertex 




B+ ^ n+RO 


23.1 ± 1.0 


r,-, 7+9.2+9.0 
'-6.0-6.9 


17.0 


„. ^+13.6 (+12.9) 
^"^■3- 8.1 (- 7.8) 


17.5 


22.5 


9, R+6.1+9.8 
^^•3-3.5-2.8 


B+ ■K°K+ 


12.9 ±0.6 


19 c;+4.7+4.9 
^^■'-'-3.0-3.8 


10.2 


,„ „+10.0{+ 7.0) 
-L-^-a. 5.6 (- 4.2) 


10.1 


12.8 


iq 1+3.7+6.0 
lo.l_2 i_j 3 


B° TT-K+ 


19.4 ± 0.6 


,Q o+T.9+8.2 
J^»-'J-4.8-6.2 


14.2 


on q+15.6(+11.0) 
8.3 (- 6.5) 


14.8 


19.2 


90 c;+5. 2+10.4 
^"•3-3.0-3.0 


BO ^ „OkO 


9.8 ± 0.6 


o 0+3.8+3.8 
^■°-2.2-1.4 


5.7 


„ ,+ 5.6(+ 5.1) 
" - 3.3 (- 2.9) 


6.3 


8.3 


o Q+2. 1+5.0 
"•^-1.3-1.1 


B° 7r-7r+ 


5.16 ±0.22 


7 n.+0.4+0.7 
' ■'^-0.7-0.7 


7.0 


P.+ 6.7 (+ 2.7) 
"■3- 3.8 (- 1.8) 


6.7 


7.1 


g g+3. 3+1.1 
"•"-1.3-0.3 


B+ -S> TT + TfO 


5.59 ±0.40 


^ q+-2+l.4 

3-»-l.l-.l.l 


3.5 


4 0+ (+ 1-7) 
1.9 (- 1.2) 


4.3 


4.1 


p. C-+2.3+1.3 
3-3_i.i_o.4 


BO ^ TT^TrO 


1.55 ±0.19 


+ 1.0+0.7 
-^■-^-0.4-0.3 


0.12 


q „„+0.50(+0.13) 
"■^^-0.20 (-0.08) 


0.1 


0.3 


1 r,+0.3+0.3 
^•"-0.1-0.1 


Acp{n+KO) 


0.009 ± 0.025 


n nn9s+O 0003+O.0009 

U.UUZO_Q QQQ3_g QQJQ 


-0.01 


-0.01 ± 0.00 (±0.00) 


-0.008 


-0.006 


n nn7+0-002+0.003 

U.UU(_Q QQJ_Q QJ^3 


AcpiTT°K+) 


0.050 ± 0.025 


n n/IQ+0-039+0 044 
U.UtM_Q Q2i_o.o54 


-0.08 


„ „,+0.03 (+0.03) 
"■"^-0.05 (-0.05) 


-0.107 


-0.053 


n nis+0.0l4+0.022 
U.UiO_Q 004-0.020 


Acp{7T-K + ) 


-0.098 ±0.012 


n n74+0-01''+0.043 
^•"'^-0. 015-0. 048 


-0.12 


n nq+0-06 (+0.04) 

"■"^-0.08 (-0.06) 


-0.131 


-0.118 


r, 1 oi +0.009+0.022 
U.iOi_0 003-0.004 




-0.01 ± 0.10 


n infi+0-027+0.056 

U.lUU_g Q3g_Q 


-0.02 


„ ^17+0. 03 (+0.01) 
"■"'-0.03 (-0.01) 


-0.001 


-0.052 


n 1 QC+0. 003+0. 004 
u.iOO_o 006-0.007 




0.58 ± 0.17 




0.70 


„ 7„+0.03 (+0.01) 
"■ '3-0.02 (-0.01) 


0.703 


0.699 


n 7t;f!+0. 002+0. 002 
"• '33-0.004-0.007 


Acp(TT~n+) 


0.38 ± 0.06 


n 1 70+0-013+0.043 
U.l(U_g QJ2-0.087 


0.14 


„ ,„+0.20 (+0.07) 
"■ °-0.12 (-0.06) 


0.206 


0.225 


n 9(50+0.043+0.059 
U.ZDU_o 032-0.063 


Acp{7T+7t'') 


0.06 ± 0.05 


-0.0002 


0.00 


0.00 ±0.00 (±0.00) 


0.000 


-0.000 


-cooilfllg 




0.43 ± 0.25 


n t;7o+0-148+0.303 
U.O(Z_g 208-0.346 


-0.04 


„ ^„+0.35 (+0.09) 
"•"3_o.34(-0.15) 


-0.600 


0.382 


n /loc+0. 061+0. 169 
U.*OO_0. 025-0. 070 




-0.61 ± 0.08 




-0.34 


„ .„ + 1.00 (+0.05) 
"■^3_o.56 (-0.05) 


-0.467 


-0.504 


r, C-9.+0. 017+0. 003 
3-3^^-0. 004-0. 017 



6.2 B ^ PV decays 

Power corrections to Oj for B — t- PV, VV are not the same as that for B — t- PP as described 
by Eq. ()14p . From Tables [8lll3t we can see that an enhancement of aj is needed to improve 
the rates of S — )• pK, loK, However, it is constrained by the measured rates of p'^n and 
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Table 7: The same as Table [B] but for B KK decay modes. 



Mode 


Datal3T] 


QCDF[n] 


pQCD 


This work 


LO 


LO NLO+Vertcx NLO(a'=-f-f , 6»") 


B+ -s> K+K'^ 
BO ^ rORO 
B° K+K- 


1.36 ± 0.28 
0.96 ± 0.20 
0.15 ±0.10 


1 O+0.9+0.7 
^■"-0.4-0.3 
9 1+1.0+0.8 
^■-^-0.5-0.5 
r, in+0.03+0.03 
J^^-0.02-0.02 


1.65 
1.75 


1 o -17 1 7+0.5+0.4 
l.O l.( -^-'-0.2-0.1 

l.z 1.0 1-O_o.3_0.2 

0.08 0.09 0.09lo:Jlo.l 


Acp{K+K°) 
Acp(KOkO) 
Acp{K+K-) 


0.12 ±0.17 
-0.58 ±0.7 


n nfiA+O'008+0.018 

U.UOt_g QQg_Q QJ^g 

n inn+O-OO'i'+O-Oio 

'-'■^'^"-0.007-0.019 




130 101 im+OMT+omi 
U.IOU U.lUl 'J-J^'JJ^-0. 013-0. 069 

0.000 0.000 0.000+0+0 
-0.182 -0.184 -0.184j:H^^ti;:"o?2 



Table 8: The branching ratios (in units of 10 ^) and direct CP asymmetries in penguin dominated 
B PV decays. The other captions are the same as Table HI 



Mode 


Data[3Tl 


This work 






NLO+Vertcx 




NLC^-f-fC-lO") 


NLO'=-f-f{5°) 


NLO'^^^{20°) 


B+ -> K''^-K+ 


9.9 ±0.8 


10.3 


9.0 


7.6 


9.8 


11.8 


B+ -> iC'+vr" 


6.9 ±2.3 


6.2 


5.5 


4.8 


5.9 


7.1 


5° -> K—TT+ 


8.6 ±0.9 


8.8 


8.3 


7.2 


8.9 


10.6 


BO ^ K-Ott" 


2.4 ± 0.7 


3.5 


3.6 


3.1 


3.9 


4.6 


B+ -s> <I)K+ 


8.30 ±0.65 


9.3 


6.9 


5.6 


7.6 


9.6 


B" 4>ko 


8.3 ± 1.1 


8.9 


6.6 


5.4 


7.3 


9.2 


Acp{K''°-K+) 


-0.038 ± 0.042 


-0.017 


-0.018 


-0.020 


-0.017 


-0.015 


Acp{K*+-kO) 


0.04 ±0.29 


-0.224 


-0.123 


-0.164 


-0.103 


-0.045 


Acp(K*-TT+) 


-0.23 ± 0.08 


-0.357 


-0.355 


-0.415 


-0.327 


-0.251 


AcpiK'^TrO) 


-0.15 ±0.12 


-0.067 


-0.125 


-.114 


-0.129 


-0.143 


Acp{<t>K+) 


0.23 ±0.15 


-0.022 


-0.025 


-0.028 


-0.023 


-0.020 


Acpi<l>KO) 


-0.01 ± 0.06 


















p^TT^ modes. Because of that, with the requirement of PV final states, we use smaller Wilson 
coefficient if the emission meson in factorizable emission diagram is a vector. 

The decay modes B — )■ ttK* , pK, coK, (pK are all penguin dominant channels. Just as argued 
in Ref. [33], contribution of annihilation diagram and NLO correction have been known to be 
able to remarkably improve theory prediction. The predictions in our framework also support 
such an argument, contributions from NLO are helpful for our theoretical predictions while still 
smaller than that in pQCD [33]. This may be due to the following fact: the system scale we use 
is higher than the stop scale used in pQCD [33! (when running scale is lower than this scale, 
pQCD uses all running parameters in stop scale for present scale and then this stop scale is 
more or less the typical scale of pQCD calculation); as a result the NLO correction in our scale 
is smaller than that in pQCD framework. 

For the case of i? — t- K(j) or B ^ ttK* decays, contributions of annihilation diagrams are more 
important than that in PP modes. This is due to the fact that emission penguin amplitude does 
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Table 9: The same as Table [51 



Mode 


Data 31] 


This work 






NLO+ Vertex 










B+ p+K° 


8.0 ± 1.45 




7 1 


7 1 


o.o 


u.o 


B+ ^ p°K+ 


3.81 ± 0.48 


3.0 


3.3 


2.8 


2.6 


2.4 


B" p~K+ 


8.6 ± 1.0 


5.4 


6.2 


7.3 


7.3 


7.2 


BO ^ pOj^o 


4.7 ± 0.7 


2.8 


3.9 


4.9 


5.0 


5.0 


B+ ujK+ 


6.7 ± 0.5 


2.4 


3.6 


4.3 


5.3 


5.2 


B° LdRO 


5.0 ±0.6 


1.9 


3.2 


4.1 


4.8 


4.6 


Acp{p+K°) 


-0.12 ±0.17 


0.016 


0.013 


0.014 


0.014 


0.014 


Acp{p°K+) 


0.37 ± 0.11 


0.635 


0.727 


0.594 


0.463 


0.285 


Acp{p-K+) 


0.15 ±0.06 


0.605 


0.549 


0.373 


0.290 


0.196 




0.06 ±0.20 


0.056 


-0.136 


-0.044 


-0.015 


0.015 


Acp{ujK+) 


0.02 ± 0.05 


0.453 


0.404 


0.167 


0.091 


0.015 


Acp{ujK") 


0.32 ±0.17 


-0.011 


0.117 


0.048 


0.026 


0.002 



not receive contribution of ag but only from 04. Our predictions without annihilation strong 
phase seem to approach the data, while the CP violation is not so good. Therefore, a small 
strong phase 9*^ = 5° is enough just like in ttK decays. 

Predictions of i? — t- pK decay modes are smaller than experimental data, because emission 
amplitudes of these channels are attributed to the destructive combination of the Wilson co- 
efficients 04 — 2^x0,%-, and hence are almost canceled out; the decay amplitudes are therefore 
dominated by annihilation contributions. Adding a strong phase 9°" = 60° for annihilation am- 
plitude will remarkably improve the branching ratio, especially the CP asymmetry. Moreover, 
the difference between Br(/9^i^*') and 'Bt{p^ K^) is accommodated with this effect. 

The case of i? — t- ujK is similar to B ^ pK, where emission diagram is also almost canceled 
out. In pQCD framework [33], vertex correction to 05 might enhance branching ratio of this 
channel in lower scale diagram Pc (color-suppressed penguin diagram), but this enhancement 
is not obviously in our framework since we use uniform scale for all diagrams. If considering 
correction of effective Wilson coefficient, we could note that the effect of larger 05 is obvious. 
Just as in S — )• pK decay mode, we add a strong phase for annihilation amplitude to improve 
predictions of branching ratios and CP violations. 

6.3 B ^VV Decays 

Naive factorization without annihilation contribution predicts a longitudinal polarization frac- 
tion near 100% for all B — )• VV decay modes, while the polarization anomaly (/l is about 50%) 
is observed by the BarBar [35) . Belle [36] and CDF \37\ experiments. Furthermore, experimen- 
tal data on B ^ pp channels show that Br(p+/9^)~ 3v{p'^p^) which breaks the factorization 
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Table 10: The same as Table |8] but for tree dominant B — > PV decay modes. 



Mode 








This 






























dcfciult 


(■—45° n°\ 
\ 'iu ,u ; 


O") 


(n° —45°') 


(n° 45°') 


B+ p+tt" 


1 n Q -1- 1 ^ 


12 


1 Q 


14 2 


lO.O 


1 7 


14 2 


B+ 


Q O _l_ 1 O 

o.o in j-.o 




7 A 


7 


7 >^ 


7 9 


7 4 




15.7 i 1.8 


19.6 


17 A 


16.5 


19.1 


17.5 


17.4 




7.3 lb 1.2 


6.2 


6.5 


6.5 


6.5 


6.1 


7.5 


B° p°n° 


2.0 ± 0.5 


0.2 


1.3 


1.5 


1.2 


1.5 


1.1 


B+ ^ K-OK+ 


0.68 ±0.19 


0.6 


0.3 


0.3 


0.3 


0.3 


0.3 


BO ^ K-ORO 


< 1.9 


0.6 


0.4 


0.2 


0.8 


0.8 


0.1 


Acpip+TT") 


0.02 ±0.11 


0.255 


0.199 


0.196 


0.133 


0.195 


0.131 


Acpip°n+) 


o.istn? 


-0.308 


-0.344 


-0.330 


-0.269 


-0.309 


-0.285 


A-Cp(p'*'Tr~) 


0.11 ± 0.06 


0.120 


0.126 


0.108 


0.066 


0.121 


0.127 


Acp(p~'n-+) 


-0.18 ± 0.12 


-0.281 


-0.282 


-0.283 


-0.281 


-0.217 


-0.176 


Acp{p°nO) 


-0.30 ± 0.38 


0.058 


0.187 


0.112 


0.381 


0.258 


-0.008 


AcpiK*OK+) 




0.191 


0.257 


-0.342 


0.205 


0.112 


0.837 


Acpik'-OR'^) 




0.000 


















Table 11: The branching ratios (in units of 10^®) and CP asymmetries for tree dominant B — > Pydecays. 
Predictions of pQCD method correspond to NLO and vertex correction results (LO results are listed in 
parenthesis). Here we use strong phase 9'^p—5° and 9^^=60°. 



Mode 


Datai31. 




pQCD(LO)|33l|34] 


this work 
LO NLO+Vetex NLO(a<=-^-^) 


B+ ^ p+ttO 
B+ ^ pOn+ 
B" -s> p+TT- 
B° -s> p-7T+ 
fiO ^ pOnO 


10.9 ± 1.5 
8.3± 1.3 
15. 7± 1.8 
7.3± 1.2 
2.0 ± 0.5 


11 0+1.8+1.4 

7+2.7+1.7 
».'_1.3_.14 

It; Q+1.1+0.9 
10-y_l. 5-1.1 

Q 9+0.4+0.5 
y-^-0.7-0.7 

1 0+1.7+1.2 
^"^-0.6-0.6 


6 ~ 9 
5 ~ 6 

0.07 ~ 0.11 


12.0 12.0 13.9+32+09 
5.4 5.2 7.4l?:^lJ;02 
18.6 19.6 17.4+42+09 
6.9 6.2 6.5l|gl°;; 

0.2 0.2 i-siElitSI 


B+ -5- K*°K+ 
-> K*0K° 


0.68 ± 0.19 
< 1.9 


r, on+0.20+0.31 
^■"'^-0.17-0.28 
r, 47+0.36+0.43 
-0.17-0.27 


S9+0-12 

"■■'^-0.07 
49+0.15 
^^■^^-0.09 


0.3 0.6 0.3lHl"o:^ 
0.3 0.6 0.4t[;;Jtg;2 


Acpip+TrO) 
Acp{p°7r+) 

Acpip'^TT-) 
Acpip^TT^) 

Acpip"^") 


0.02 ± 0.11 

o.i8«:?? 

0.11 ± 0.06 
-0.18 ±0.12 
-0.30 ±0.38 


n 007+0 021+0.080 
'J-^y -0.031-0.103 
n 008+0 034+0.114 
'J-^'aO-0. 026-0. 104 
n nzlA+0-003+0.058 

U.Uiy:_Q QQ3_Q ggg 

r, 997+0.009+0.082 
'J-^^' -0.011-0.044 
r, -npi+0. 050+0. 235 
U.ilU_g Qjy_Q 


0-20 
-20 ~ 

-75 ~ 


0.251 0.255 0.1991^4110.053 
-0.351 -0.308 -0.344+°°f +°0°g^ 

0.113 0.120 o.uet^^'Z 

-0.225 -0.281 -0.284+g-064+0.045 
0.048 0.058 0.187lg:»tl°or7 


Acp{K*Ok+) 
Acpik*OKO) 




n nsQ+0 011+0.028 

U.UOM_g Qjj_Q 

n n'?t;+o oi3+o.oo7 

U.UOO_Q Q]^7_Q 020 


n nf;Q+0.056+0. 010+0. 092+0. 040 

U.UDM_0 Qg3_Q 003_0 Qgg_Q 00gQ 


0.360 0.191 0.257lg:°42lo.o?| 
0.000 0.000 oIqIq 
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Table 12: The branching ratios(in units of 10^^) and CP asymmetries for penguin dominant B — > 
PT^decays. Predictions of pQCD method correspond to NLO and vertex correction resuhs (LO resuhs 
are hsted in parenthesis ). Here we use strong phase 0^j^t=0'^,=5°. 



Mode 


Datalsn 


QCDF[TT1 


pQCD(LO) 33 


this work 


LO NLO+Vcrtcx NLOia^ff ,6'^) 


B+ ^ K''0^+ 
B+ ^ K'+wO 
gO ^ K'-n+ 


9.9 ± 0.8 
6.9 ± 2.3 
8.6 ± 0.9 
2.4 ± 0.7 


+0.7+2.4 
' -0.7-2.2 
Q 9+1.0+3.7 
^■^-1.0-3.3 

o c;+0.4+1.6 

■J-3_0.4-1.4 


6.0t?;«l?;I(5.5) 

4.3t5;»2ll-i(4.0) 
6.0t|Jl?-j(5.1) 
2.0tJ;n;l(l-5) 


7.5 10.3 9.8t.f,rlLl'g 
4.7 6.2 5.9l?°;j;^ 
6.5 8.8 8.9l?:ll?* 

9 c o c: r> n + 0. 8+1.1 
z.o o.o •J-»_o.6-0.3 


B+ <t>K+ 


8.30 ± 0.65 
8.3± 1.1 


a 0+2.8+4.4 
o-o_2. 7-3.6 
o 1+2.6+4.4 
"■^-2.5-3.3 


7.8t5-|+5-8(i3,8) 
7-3t?11'ii(12.9) 


10.8 9.3 7.6l}:|lg| 
10.4 8.9 T-^tlitoi 


Acp{K*07T+) 

Acp{K' + 7T-) 

AcpiK'OnO) 


-0.038 ± 0.042 
0.04 ± 0.29 
-0.23 ±0.08 
-0.15 ±0.12 


n nnA+'' 0i3+0.043 

'JlJU't-0. 016-0. 039 

n nifi+0-03i+o.iii 

'-'■'^^"-0.017-0.144 
p. 1 91 +0.005+0.126 
U-1^1_0. 005-0. 160 

n ins+o ois+o.09i 

u.luo_Q 028-0.063 


-o.oiiroSl'^:SJ(-o-o3) 

-0.32ir2^t°o:f9(-0.38) 
-0.60l°:f,lg:?°(-0.56) 
-0.1liHIirol(-0.60) 


-0 021 —0 017 Q Qi -7+0.003+0.012 
U.UZi U.UK 'J-'J^ ' -0.002-0.003 

-0.348 -0.224 -O.lO3loo27lo.oi6 

-0 443 -0 357 -0 327+'' °^°+"-°^° 
U.00( '-'■•J^ ' -0.008-0.016 

0.012 -0.067 -0.130in28to.M9 


Acpi^K+) 
Acp{4>K°) 


0.23 ± 0.15 
-0.01 ± 0.06 


n nnfi+'^''00i+0'00i 

U.UUU_Q QQ1_Q 


o.oit?;.!]?tro?(-o-02) 

0.03l-°^l°:?°(0.00) 


-0.022 -0.022 -O.O23ino2lo.m4 

0+"+° 
u u "-0-0 



Table 13: The same as Tabled but with ^°a-=^pA'=60°- 



Mode 


Data'ST 


QCDF[TT] 


pQCD{LO) 33 


this work 


LO NLO+Vertex NLO(a<=-''-'', 6*") 


B+ -s> p+K° 
B+ -s> p''K+ 
-s> p-K+ 
BQ _+ pO^o 


8.0 ± 1.45 
3.81 ± 0.48 

8.6 ± 1.0 

4.7 ± 0.7 


7 0+6.3+7.3 
' ■°-2.9-4.4 
o c;+2.9+2.9 
■^■3-1.2-1.8 
Q R+5.7+7.4 
o-D_2.8-4.5 
t- 4+3.4+4.3 
°-^-1.7-2.8 


8.7ll«ir3(3.6) 

5-ll2.8-2.6(2-"') 

8.81111^3(4.7) 
4.8l^;|l|:g(2.5) 


4.2 5.2 e.Sloali'? 

2.3 3.0 2.6i;^ii5;^ 

4.9 5.4 7.3+°qI+11 
2.7 2.8 5.0ir4li:^3 


B+ ujK+ 
BO ^ ujRO 


6.7 ± 0.5 
5.0 ±0.6 


4 0+4.4+3.5 
^■"-1.9-2.3 
4 1+4.2+3.3 
^■1^-1. 7-2. 2 


10.6lJV-y4'{2.1) 
9.8im.I(l-9) 


2.4 3.6 5.31^41^:^5 
1.9 3.2 4.8l»-l- 


Acp{p+KO) 
Acp{p''K+) 
Acp(p-K+) 
Ac pip'' K^) 


-0.12 ±0.17 
0.37 ± 0.11 
0.15 ±0.06 
0.06 ±0.20 


n nn'^+0-002+0.005 

U.UUO_Q QQ3_Q 002 

n 4t:4+0. 178+0. 314 
"•^'^^-0. 194-0. 232 

n qiQ+0. 115 + 0. 196 

U.ol3_o. 110-0. 127 
n nH7+0-012+0.087 
u.uo/ _Q oi2_o.o68 


o.oi1H;1°o:°o1(o.o2) 
o.7ilg:il°o:lI(o.79) 

0.64l»:^^lHI(0.83) 
0.07l°:rroI(0.07) 


0.019 0.016 0.014l°;»;i»o.Sl^ 
0.726 0.635 0.463lo.o36lo.oi4 
0.593 0.605 0.290l°;°^J120i2 
-0.040 0.056 -0.015l°;°g510 025 


Acp{ujK+) 
AcpiujK°) 


0.02 ±0.05 
0.32 ±0.17 


n 991 +0.137+0.140 
U.ZZ±_o i2g_o 130 

n nA7+0-01*+'3-055 
U.Ut(_Q 016-0.058 


0.32l»;i|l»o:°ot(0.32) 
-0.03lH^iro3(-0.03) 


0.688 0.453 ■09ltl°ll1:°o°o7l 
0.065 -0.011 0.026lg;°°ll°;°°| 
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Table 14: Branching ratios for B VV decay modes (in unit of 10 ^) which includes the contribution of 
effective Wilson coefficients and effect of different strong phase 0° = 60° ± 15°) for annihilation diagram. 





Mode 


Data jS^l^ 


This work 








NLO+Vertex 




NLO'=^'^(45°) 


NLO=-''-''(60°) 


NLO'= (75°) 


B+ 




24.0 ± 2.0 


13.4 


16.8 


16.8 


16.8 


16.8 


B° 




24.2 ± 3.1 


22.3 


19.8 


21.7 


22.3 


22.7 




-> pOpO 


0.73 ±0.27 


0.4 


0.92 


0.67 


0.61 


0.57 


B+ 


K*Op+ 


9.2 ± 1.5 


16.2 


14.0 


9.6 


8.3 


7.2 


B+ 


-> K *+pO 


< 6.1 


9.9 


9.0 


6.4 


5.6 


5.0 


BO 


_^ ^•+p- 




13.9 


13.0 


9.1 


7.9 


6.9 


BO 


^ A-'OpO 


3.4 ± 1.0 


5.6 


5.2 


3.6 


3.1 


2.7 


B+ 


-s> K*°K''+ 


1.2 ±0.5 


0.9 


0.8 


0.6 


0.5 


0.4 


B° 




1.28 ± 0.35 


0.8 


0.6 


0.5 


0.5 


0.5 


BO 


-s> K'+K*- 


< 2 


0.07 


0.07 


0.07 


0.007 


0.07 


B+ 


<f>K*+ 


10. 0± 1.1 


19.4 


15.2 


10.9 


9.5 


8.4 


BO 


(pK'O 


9.8 ±0.7 


18.7 


14.8 


10.5 


9.2 


8.1 


B+ 


-> ujK*+ 


< 7.4 


5.6 


4.2 


3.3 


3.0 


2.8 


BO 


-> ujK*o 


2.0 ±0.5 


6.2 


4.1 


2.8 


2.5 


2.2 



power-counting while the power-counting Br{p^ p^)^3r(p^ p^) still satisfies. Motivated by the 
anomaly and the unsatisfied power counting relations, we shall study in detail the polarization 
in i? — 7- VV decays, especially in i? — t- pK* and (pK* decays, and the branching ratios of B ^ pp 
modes in this section. The relevant predictions are listed in Tables [T^19[ 

Before moving ahead, one point that should be noted is that we use a bigger gluon infrared 
cut-off in annihilation diagram, Jig =0.52GeV ~ V^Pg, which is constrained by the huge penguin 
diagram contributions. 

It is realized that the trick used in vrvr puzzles is no longer useful here. If we add effective 
(^1,2 just as we do in vrvr mode, it is only helpful for branching ratios of p^p^ and p^p^ modes, 
but at the cost of overflowing that of p^p^. Branching ratios of — J- p^p^ can be decreased by 
introducing proper strong phase (such as 9"" = 60°) in annihilation amplitudes, while be increased 
by considering effective Wilson coefficients a^^ 2 • Moreover, these two elements make predictions 
of branching ratios of B ^ pp, pK* , (pK* and ojK* to be better consistent with the experimental 
data. The effect of different strong phases to branching ratio, longitudinal polarization and 
direct CP asymmetry are listed in Tables [^11161 As for the channel B — )• K*K*, predictions of 
the branching ratios are smaller than the experimental data. The strong phase effect does not 
increase the predictions, so the prediction without strong phase are better than the one with it. 

Another important point should be noted that the predictions for the branching ratios of 
B — )• pK*, (j)K* and uK* modes are all bigger than the experimental data, which means that 
if we want to solve the observed polarization anomaly, we need to find some way to reduce the 
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Table 15: The same as Table [Til but for the longitudinal polarization fraction. 



Mode 


Data[3Tl 


This work 






lNijiJ+ Vertex 


IN ijij 


MT f f f Ar^O\ 

IN ijij ) 


IN i^yj (DO j 


MT f^G f f { '7 f^.O \ 
INLU ^ ^ \ (d ) 


/z,(p+pO) 


0.950 ± 0.016 


0.94 


0.92 


0.95 


0.95 


0.95 


fr 

J iv \r H ) 


0.978 ± 0.023 


0.95 


0.95 


0.95 


0.95 


0.95 


J J-i \H H J 


0.75 ± 0.15 


0.84 


0.86 


0.77 


0.74 


0.71 


fr (K*°p+) 


0.48 ± 0.08 


0.85 


0.82 


0.57 


0.45 


0.32 


fr (K*+p'>) 


0.96+J!?« 


0.86 


0.85 


0.65 


0.56 


0.47 


fL{K*-p+) 




0.81 


0.80 


0.57 


0.46 


0.34 




0.57 ± 0.12 


0.78 


0.75 


0.48 


0.36 


0.22 




75+016 
U-'0_o.26 


0.85 


0.81 


0.60 


0.49 


0.37 




0.80 ±0.13 


0.83 


0.63 


0.60 


0.53 


0.46 


fL(K*+K*-) 




0.99 


0.99 


0.99 


0.99 


0.99 


fL{4>K* + ) 


0.50 ±0.05 


0.87 


0.83 


0.58 


0.45 


0.31 




0.480 ± 0.030 


0.87 


0.83 


0.58 


0.45 


0.31 




0.41 ± 0.19 


0.90 


0.86 


0.68 


0.58 


0.48 




0.70 ±0.13 


0.93 


0.89 


0.67 


0.53 


0.37 



longitudinal amplitude and enhance the transverse ones simultaneously. Many studies have been 
made to try to provide possible resolutions to the anomaly both within the SM [351 133 SOI SI] 
and in various new physics models [12 HSl |311 US] . 

In our framework, sizable annihilation contributions appear in all the three parts of the 
amplitudes: longitudinal, parallel and transverse. In order to get a reconciled prediction for the 
polarization part, we shall add an extra strong phase {9^ = 60°) in all annihilation diagrams. As 
a result, the longitudinal amplitude can be partly counteracted by annihilation contributions, 
while the parallel and transverse ones are not be influenced much. From Table [T^ we can see that 
our predictions for the longitudinal polarization fraction are generally consistent with the current 
data when considering complex annihilation contributions and effective Wilson coefficients. 

In i? — )• VV decays, there are another two interesting observables (f)\\ and c/)^, which are 
defined to be the relative phase between the parallel and the longitudinal amplitudes and that 
between the transverse and the longitudinal amplitudes, respectively. The corresponding defi- 
nition can be found in Ref. j^41j. To compare with the data, we only discuss B — t- (pK*^ , (I)K*^ 
decay modes. Since (/)|| is the same as (p±, we simply write them as (p and the numerical results 
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Table 16: The same as Table but for direct CP asymmetries. 



Mode 


DatafST] 


This work 






ATT 1 T r J 

NLO+Vertex 


NLO'^-'^-'^ 


ATT fff/^r'0^ 

NLO ■' (45 ) 


NLO (60 ) 


ATT /~\ P f t /'-Tr'OX 

NLO^ ^ (75 ) 




—0.051 lb 0.054 


0.001 


0.001 


0.001 


0.001 


0.001 


A.n p io^ 0~ ) 


0.06 ± 0.13 


—0.002 


-0.029 


-0.041 


-0.038 


-0.033 






0.702 


0.177 


0.350 


0.417 


0.475 




-0.01 ± 0.16 


—0.005 


-0.006 


-0.009 


-0.009 


-0.008 




w.^u q 29 


0.184 


0.182 


0.266 


0.273 


0.257 


ArpiK*- n-^\ 




0.148 


0.151 


0.231 


0.231 


0.206 


AcpiK*OpO) 


0.09 ±0.19 


-0.090 


-0.101 


-0.147 


-0.155 


-0.154 


Acp{K*OK'+) 




0.081 


0.085 


0.141 


0.143 


0.128 


Acp{K*°K*°) 



















Acp{K*+K—) 




-0.261 


-0.261 


-0.261 


-0.261 


-0.261 


Acp{<PK*+) 


-0.01 ± 0.08 


-0.003 


-0.003 


-0.0006 


-0.007 


-0.007 


Acp{4>K*°) 


0.01 ± 0.05 

















Acp{ujK''+) 


0.29 ±0.35 


0.341 


0.383 


0.534 


0.522 


0.463 


Acp{ujK*0) 


0.45 ±0.25 


0.078 


0.116 


0.170 


0.182 


0.185 



are as follows: 



NLO+Vertex : 




) = 99.0°, 




= 100.7°, 


with 6" : 




) = 63.7°, 


<A(0K*°) 


= 66.1°, 


with e^and a^ff : 




) = 73.1°, 




= 75.9°, 


QCDF[llj : 




= (8oii)°, 




= {78titr, 


expt : 




= (46± 10)°, 




= (44t?)°, 






= (40± 10)°, 




) = (43 ± 7)° 



It is noted that, although our predictions are bigger than the data, they are roughly consistent 
with the results in QCDF 06]. 



7 Conclusions 

Based on the approximate six-quark operator effective Hamiltonian derived from perturbative 
QCD, the QCD factorization approach has been naturally applied to evaluate the hadronic 
matrix elements for charmless two body B-meson decays. It is shown that, with annihilation 
contribution and extra strong phase, our framework provides a simple way to evaluate the 
hadronic matrix elements of two body decays. 

For B —7- PP final states, our predictions for branching ratios and CP asymmetries are 
generally consistent with the current experimental data within their respective uncertainties. 
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Table 17: Branching ratios for B — > VV decay modes (in unit of 10^^). The central values are obtained 
with /ig=0.52GeV and 6°- = 60°. The first error arises from the varying for fiscaie — 1-4 ^ 1.6 GeV, the 
second one stems from the shape parameters of light mesons. 



iViOQt; 








this work 










T O 






L\Li\^y(l ^ ^ ^ u ) 




9zL n -1- 9 n 


9n n+4-0+2.0 


17+9+ 1 
1 ( in z m 1 


1 Q V 


1 Q /I 

10.4 


1 Q /I 
10.4 


IK 0+9.9+1.2 
lO-S-4.4-0.7 


n —^pp 


24 2 it 3 1 


c + 1.5+2.4 
zO.O_2 g_i 5 


oo in u in ^ 


21.1 


22.3 


24.9 


gr, ri + 14.0+1.3 

^^■^-6.2-0.7 




0.73 it 0.27 


"■^-0.4-0.2 


0.9 it 0.1 it 0.1 


0.3 


0.4 


0.4 


n fj+0. 2+0.1 
'-'■"-0.1-0.0 


B+ -s> K*°p+ 


9.2 it 1.5 


Q 9+1.2+3.6 
="■^-1.1-5.4 


17 (13) 


11.4 


16.2 


9.8 


o 0+3.4+3.5 
"■°-1.2-2.4 


B+ -s> K*+p° 


< 6.1 


c P.+0.6+1.3 
■^■"^-0.5-2.5 


9.0 (6.4) 


7.3 


9.9 


6.2 


c- Q+2. 1 + 1.9 
^■^-1.0-1.2 


B" ^ K*+p~ 


< 12 


o „+l.l+4.8 
»-y_1.0-5.5 


13 (9.8) 


10.2 


13.9 


8.5 


o 0+2.1+3.2 
"■•^-1.0-2.5 


BO ^ K*0pO 


3.4 it 1.0 


A fi+0.6+3.5 
4-0_o.5-3.5 


5.9 (4.7) 


3.9 


5.6 


3.4 


q O + 0.5+1.7 
■^■•J-0. 2-1.1 


B+ -5. 


1.2 it 0.5 


n C.+0.1+0.3 
'^■"-0.1-0.3 


0.48 


0.7 


0.9 


0.6 


n c:+0.2+0.2 
^■^-0.1-0.1 


B° -5- K*°K''° 


1.28 it 0.35 


n 0+0.1+0.2 

'-'■"-0.1-0.3 


0.35 


0.5 


0.8 


0.6 


n c:+0.2+0.2 
'^■^-0. 1-0.1 


_^ K*+K*- 


< 2 




r, 1+0.0 + 0.1 

'^■-'--0.0-0.1 


0.07 


0.07 


0.07 


r, ny+O.Ol + O.OO 

'^■"' -0.01-0.01 


B+ -5- !j>K*+ 


10.0 it 1.1 


^'-'■^-1.3-6.1 


15.96 


15.9 


19.4 


12.4 


Q ^+2.5+2.4 
^■"-0.6-1.6 


B° -5- !j>K*° 


9.8 it 0.7 


Q ^+1.3+11.9 
^■'-'-1.2-5.9 


14.86(10.2^2;^) 


15.4 


18.7 


11.8 


Q 9+2.3+2.3 
^■^-0.5-1.6 


B+ LuK*+ 


< 7.4 


.1 n+0-4+2.5 
'^■'-'-0.3-1.5 


7.9(5.5) 


5.4 


5.6 


3.7 


q n+l. 0+2.1 
■J-U_o. 4-1.0 


B° -> LuK*° 


2.0 it 0.5 


r, c:+0.4+2.5 
^■■^-0.4-1.5 


9.6(6.6) 


5.8 


6.2 


3.8 


r, c: + 0. 7+1.0 



Table 18: The same as Table [TTl but for longitudinal polarization fractions. 



Mode 


Data[3T] 


QCDF liij 


pQCD Ha [48] 


this work 










LO 


NLO+Vertex 


NLO(6l") 


NLO(a'=^^,e'') 


fL{p+P°) 


0.950 it 0.016 


n 013+0.01+0.02 

'^■™-0.01-0.02 


0.94 


0.95 


0.94 


0.95 


n Qcr+o.oo+o.oo 

'-'•='"'-0.00-0.00 


fL{p+p-) 


0.978 it 0.023 


r, Q9 + 0.01 + 0.01 

'^■^'^-0. 02-0. 02 


0.94 


0.96 


0.95 


0.95 


n Qc+o.oo+o.oo 

'-'•='"'-0.00-0.00 


fL{p"d') 


0.75 it 0.15 


r, 00+0.03+0.06 
^^■^^-0.04-0.37 


0.60 


0.79 


0.84 


0.56 


n 7-+0.03+0.01 
'-'• '^-0.03-0.01 


fL{K*Op+) 


0.48 it 0.08 


n 40+0.03+0.52 
"J-^°-0.04-0.40 


0.82 (0.76) 


0.79 


0.85 


0.51 


n .c+O.05+0.02 
'-'•^"'-0. 01-0. 06 


fL{K* + p°) 




'J-yO-o.io 


n C7+0.02+0.31 
'^■'-" -0.03-0.48 


0.85 (0.78) 


0.83 


0.86 


0.59 


n c-7+0.05+0.02 
'-'•■"-0.04-0.04 


fL{K*-p+) 




„ c:q+0.02+0.45 
'^■'^ •-'-0.03-0. 32 


0.78 (0.71) 


0.77 


0.81 


0.49 


n Af!+0.04+0.03 
'-'•^"-0.02-0.07 


/i(X*OpO) 


0.57 it 0.12 


r, qn+o.oo+o.eo 

'^■■'^-0. 00-0. 31 


0.74 (0.68) 


0.71 


0.78 


0.38 


n qc+0.01+0.02 

'-'•■^°-0.01-0.05 




75+0.16 
'^•"^-0.26 


r, f;r, + 0.01 + 0.42 

'-'■'^^-0.02-0.33 


0.72 


0.80 


0.85 


0.53 


n AQ+0.06+0.02 
'-'•^^-0. 03-0. 05 




0.80 it 0.13 


n cQ+0.01+0.34 
'-'■'-'^-0.01-0.27 


0.67 


0.77 


0.83 


0.57 


n cq+0.07+0.01 

'-'•"'•-'-0.03-0.03 


fUK*+K'-) 






0.99 


0.99 


0.99 


0.99 


n QQ+O.OO+O.OO 

'-'•='^-0. 00-0. 00 


M'PK*+) 


0.50 it 0.05 


n 40+0.04+0.51 
'-'■^''-0.07-0.42 


0.748 


0.80 


0.87 


0.53 


n 4cr+0.03+0.02 
'-'•^"'-0. 01-0. 04 




0.480 it 0.030 


n c:n-i-o.04+o.5i 

'-'■^"-0.06-0.43 


0.75(0.59tE;;02) 


0.80 


0.87 


0.53 


n AK+0.03+0.03 
'-'•^"'-0. 01-0. 04 




0.41 it 0.19 


(.^+0.03+0.32 
'-'■"' -0.04-0.39 


0.81(0.73) 


0.84 


0.90 


0.63 


n co+0.03+0.02 
'-'•"'°-0.01-0.02 


fL{u>K*0) 


0.70 it 0.13 


n Ko+0.07+0.43 
'-'■^"-0.10-0.14 


0.82(0.74) 


0.86 


0.93 


0.64 


cq+0.03+0.04 
"•"'•-'-0.01-0.04 
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Table 19: The same as Table [T71 but for direct CP asymmetries. 



iVloae 


Jjata oi- 






this work 










LO 


NLO 


NLO(e'') 




Acpyp^p ) 


— (J.Uol ± U.U54 


r\ r\r\r\c 

O.OUOb 





0.001 


0.001 


0.002 


U.UU±_Q 000-0.000 
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once the effective Wilson coefficients and annihilation amplitude with small strong phase = 
5°) are adopted. Especially for the branching ratio of i? — t- vr^vr'^ mode, our result, although 
being still smaller than the data, is consistent with that in QCDF [llj. 

As for B — )• PV decays, similar conclusions are found. The exceptions here are the branching 
ratio of p+vr'^ and B — t- K*^7r^ modes, which are bigger than the data. An interesting point 
should be noted that our predictions are also consistent with the ones in QCDF [TlJ. Since 
the current data on CP asymmetries have large uncertainties in these modes, more precise 
experimental data are expected to further test our framework. 

In 5 —7- VV decay modes, we have shed light on the polarization anomalies observed in S — ^ 
(j)K* , pK* ,u}K* decays. It is noted that these anomalies could be explained in our framework 
when considering annihilation contributions with a strong phase {O"" = 60°). Moreover, the 
annihilation contributions with a strong phase have remarkable effects on the branching ratios 
and CP asymmetries, especially on the observables of penguin dominated decay modes. As a 
result, with the effects, we also have good predictions for the branching ratios and direct CP 
asymmetries. 

It is noted that the method developed in this paper allows us to calculate the relevant 
transition form factors. Our predictions (for B to light mesons form factors) are consistent 
with the results of light-core QCD sum-rules and pQCD. We further apply the method to 
Bg — )• PP, PV, VV decays, and the paper |50J is now in preparation. 
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Appendix: Calculations of Hadronic Matrix Elements 

With the considerations and analysis in the text, the QCD factorization approach with six- 
quark operator effective Hamiltonian enables us to evaluate all the hadronic matrix elements of 
nonleptonic two-body B-meson decays. 

When generalizing the above analysis to the present framework based on the approximate 
six-quark operator effective Hamiltonian, there are in general four types of six-quark diagrams, 
and each one corresponding to four types of effective six-quark operators. So their hadronic 
matrix elements for two body mesonic decays lead to sixteen kinds of diagrams (see Figs. Hl^al)- 
(d4)) as each of the effective six-quark operators leads to four kinds of amplitudes in the QCD 
factorization approach. In order to get a global form of hadronic matrix elements, we should 




Figure 4: Four types of effective six quark diagrams lead to sixteen diagrams for hadronic two 
body decays of heavy meson via QCD factorization. 

rewrite mesons with SU(3) symmetry into a vector/matric form. We use a column vector Ap 
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containing CKM matrix elements, collect the three B-meson states into a row vector B, and 
represent the final-state pseudoscalar and vector mesons by matrices P and V 

\ 



B 



V 



V 



V K*+ 



V2 y/2 

K+ 
K*- \ 



IT 



K- 

Vs+V's J 



K 



*o 





( ° 
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Ap = 













(A-1) 



With the above definition, we can express hadronic matrix elements with the sixteen diagrams 
in FigH] marked with (al)-(d4): 

A{B M1M2) = 

^ ^ [Ci{Mf + M^^ + Mf + Mf^)BMrUiM^Kp 

2=1,4,6,8,10 p=M,c 

+ Ci{Mf + Mf + Mf + Mf^)BMiAp ■ It [C/iMa] 

+ Q(Mf3 + Mf + Mf + Mf)BUiMxM2Ap 

+ Ci{Mf + Mf + Mf + Mf^)BAp ■ Tr [C/iMiMa] | + 

Yl {Ci{M^^ + Ml'^ +Mf^ + Mf^)BMiAp-TT[UiM2] 

j=2,3,5,7,9p=«,c 

+ Ci{Mf + Mf + Mf + Mf)BMxUiM2 Ap 
+ Ci{Mf + Mf + Mf + Mf)BAp ■ Tr [[Z^MiMa] 
+ Q(Mf + Mf + Mf + Mf )SC/,MiM2Ap 

with the matrices C/,- defined as 



(A-2) 



/ 5p« \ 


\ y 



J7p fori = 1,2, 
/ for i = 3 — 6 , 
Q for i = 7 - 10 

/ 1 \ 

10, Q 

Vo 1/ 







' 3 





\ 



\ 



4/ 



(A-3) 



Definitions of Mj are attributed to three types of current-current four-quark vertexes: 



M 



Mll for i = 1 -4,9,10 {V -A)x{V -A), 
Mlr fori = 5, 7 {V - A)x{V + A), 

MsP fori = 6, 8 {S-P)x{S + P). 



(A-4) 
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Therefore, there are totahy 48 kinds of hadronic matrix elements involved in the QCD factoriza- 
tion approach, while it is easy to check that only half of them are independent with the following 
relations: 
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(A-5) 



where T^^ and T^j^ {X = LL, LR, SP) represent the factorizable emission diagram contri- 
butions, rf„ and Tf^ {X = LL, LR, SP) are the non-factorizable emission diagram contri- 
butions. A^^, A^^ and A^^, A^^ {X = LL, LR, SP) denote the so-called factorizable and 
non-factorizable annihilation diagram contributions, respectively. With this definition, we can 
rewrite Eq. ()A-2p into following form 

A[B M1M2) = 

+ B (^A^^^^^^'iB) Up + P*^i^^2(B)P + P^^'^^^B)P q) MiMaAp 
+ BAp ■ Tr \(e^'^^^'^{B) Up + P^^^^^'{Bf + P^^'^^^{B)p q) Mi Mai | ,(A-6) 



p=u,c 



BMiAp ■ Tr 



where Mi and M2 are pseudoscalar (P) or vector(y) mesons. The twelve types of ampli- 

C ^EW 



tudes r^'^iA^2(M), C^^^^^iM), P^^^\M), P^''^^^{M), A^''^^^^{M), E^'^^^{M), 
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+ [Ce 

\M) = Anas 

+ [C7 
+ [C8i 

for the so-called emission 



pMiM2i 



oCMiM2 
EW 



pMiM2 



T3AM1M2 



pEMiM2 
^EW 



(M) 



(M) 



(M) 



Ait as' 

Airas' 
+[Cq 

A-Kag- 

+ [C5 

A-KUs' 
+ [C7 
Airag' 



pM.M2^M), P^^^^M), P^^^^^(M), P|^^^^(M), Pi^^^^(M), are defined as follows 

^M^M2^M) 



A^) + ^C'r(/.)]r|p^^^^(M) + ^C,(/.)rf/^^^(M)}, (A-7) 
diagrams, and 

+ ^Crii^M^sp'^^'iM) + ^Cv(mX^^^^(M)}, (A-8) 



for the so-called annihilation diagrams. Where T^^, T^^, ^^A' i-^ ~ L,L,LR,SP, 

A = a,b) arise from the hadronic matrix elements and their detailed expressions are given 
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below. 

We now list our results for the various decay amplitudes expressed in terms of above topolog- 
ical amplitudes. The detailed calculations of the hadronic matrix elements for B — )• PP decays 
could be found in our previous paper [7J. As for decay amplitudes and the hadronic matrix 
elements for B — )• PV, VV decays, we shall list them one by one as follows. Firstly, for B ^ np 
decay channels, we have 

A(BO^vr-p+) = VM[P-p{B) + ^P^^^{B) + P^P{B) + PY{B) + P^{B) + ^pX{B) 

-V^dV:,[T^^\B) + E-^P{B)l 
A{B'^p-7r+) = Vt,V,l[PP^B) + ^P^^{B) + Pg'{B) + P^/{B) + PY{B) + lpX{B) 

-VudV:,[T'^B) + E^^B)], 
^(fi+^vrV) = -^{Vtm-P^'{B)+P'^B) + p-/^iB) + ^P^;,'{B) + ^P^{;^iB) 

_pA.p^B) + pr(^) - Pew'{b) + p^^m - v^av:,\Tp-{B) 

+CP^{B) + AP^{B) - A^'PiB)]], 
A{B+^p+^') = 1={V,,V:,[-PP^{B)+P-P{B) + PP^^{B) + \P^^^{B) + \P^^^^ 

-P^P^B) + P^iB) - P^iB) + P^^^(i?)] - V^,V:,[T-P{B) 
+C''P{B) + A^'PiB) - AP^iB)]}, 
A(fi°^7rV) = ]^{A{B^ ^ p+TT-) + A{B^ ^ p-n+) 

-V2[A{B+ ^ vr+pO) + A{B+ ^ p+vr")]}, (A-9) 

For B — )■ 'kK* decay channels, the decay amplitudes are 

A{B+ ^ 7r+K*0) = -I4,y;jP-^*(i3) - \P%'^\B) + P|^*(S) + IpE^'^^B)] + VMA-^'{B), 

a{b+^7t'k*+) = -1={vm[p-^\b) + pEw''{b) + ^pE^^*{b) + pe^\b) + ^pE^^*{b)] 

-VusV:,[T^''* (B) + C^*^(P) + 
A{B' ^ 7,-K*+) = Vt,V,l[P-^\B) + + - IpEw^'^B)] - VusV:,T-^'{B), 

A{B^^n'K^) = -1={V,,V:,[P^^\B)-PE;^{B)-^pE{^''\B) + P^^\B) 

-^P§^^\B)] + VusV:,C^'^B)}, (A-10) 
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For B pK decay channels, we have: 

A{B+^p+K^) = -VtsV,l[P^''{B)-^P^^^{B) + Pf{B) + ^P^^^{B)] + Vu,V:,A^''{B), 
A{B+^p'K+) = l^{VuV:,[P^^ {B) + pZ{B) + \P^E^ {B) + Pf {B) + '^-P^^ {B)] 

A{B' ^ p-K+) = VumiP'^^'iB) + ^P^#(5) + Pf{B) - ^pX(S)] - VMT'>^{B), 
A{B^^p'K') = -l^{VuV:,[PP^{B)-P^P,{B)-\p''^^{B) + Pf{B) 

-\PE^m + VM^^'m, (A-11) 

For B (jjK decay channels, we have: 

A{B+^uK+) = ±{VuV,l[P^''{B) + ^PE^^iB) + ^P^i^iB) + 2P^^{B) + P^^{B) 

+^P|if (B)] - VusV:,[T'^''iB) + C^-(B) + A^^{B)]}, 
AiB'^uK') = -^{VuV*,[P^^{B)-^Pg^^{B)+^-P^{^{B) + 2P^-{B) + P^^{B) 

-^Pfi^f (5)] + VusV:,C^'^{B)}, (A-12) 
For B — >■ (pK decay channels, we have: 

A{B+^cl>K+) = VuV,l[P^HB)-lpSw'^{B)-^P^i{B) + 2Pg'^{B) + Pf{B) 

- V^sV:,A^''{B), 

A{B^^^K^) = -VuVl[P''^{B)-\p^^'^{B)-\p^^{B) + 2P^^{B) + Pf{B) 

-ipX(S)], (A-13) 

For B K*K decay channels, we have: 

A{B+^K*'K+) = -VuV:,[P''^\B)-\p^^^\B) + P^^'\B) + \pE§''{B)] 

+VudV:,A^^{Bl 

A{B'^K*'K') = -VtaVl[P^'^{B)-\pq^''^{B) + Pf^{B) + Pr^{B) 

+Pr\B) - IpirHB) - lpi^^\B) - 

(A-14) 
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For B ^ pp decay channels, we have(three part of amplitude are similar): 

A{B'^p+p-) = VuV,l[P^%B) + ^P^^^{B) + Pg'iB) + 2P^^{B) + ^P^^{B) 

-VudV:,[TPP{B) + EPP{B)l 
A{B+^p+p') = -^{VuV;,[P^'^{B) + P^^^{B)]-V^,V:,[T''P{B) + C^^^^ 

A{B^^p'p°) = ^{-VtaV,l[P^P{B)-P^/^{B)-^P^^{B) + Pg'{B) + 2P^/{^^ 

+lPE^iB) - Ip^^m + V^V:,[-CPP{B) + E^^{B)]}, (A-15) 

For B — > pK* decay channels, we have(three part of amplitude are similar): 

A{B+^p+K*^) = -VtsV,l[P'^\B)-^-P^^^\B) + Pf\B) + ^PX^\B)] 

+VusV:,A'''<\B), 

A{B+^p'K*+) = ±{VuV,l[P^^\B) + P^;,^{B) + lp^^^\B) + Pf\B) 

+^pX^\b)] - VusV:,[T'>^\B) + C^>(i?) + AP^\B)]}, 

A{B^^p-K*+) = VuV:,[P^''*{B) + ^P^^^*{B) + Pf\B)-^PX^\B)] 

-VusV:,TP^\B), 

A{B'^p'K') = -l={V,,V,l[P^''\B)-P^;^'iB)-^P^^^\B) + Pf\B) 

~Pe^^\b)] + VusV:,C'''^{B)}, (A-16) 
For B — )• coK* decay channels, we have(three part of amplitude are similar): 
A{B+^coK*+) = l={VuVtl[P^''*{B) + ^PE^^\B) + ^P^^^{B) + 2Pf^{B) + P^^\B) 

+^P|if *(S)] - VusV:,[T^''\B) + C^*-(i?) + A-^*(i?)]}, 
A(B0^a;O = -l={VuVtl[P^''\B)-^PE^^\B) + ^P^;,^{B) + 2Pf^{B) + P^^*{B) 

-^P§-^\B)] + VMC'''''{B)}, (A-17) 
For B (f)K* decay channels, we have(three part of amplitude are similar): 

A{B+^<t^K*+) = \^,V^[P^''^iB)-^P^^'^{B)-^P^;^iB) + 2Pf^{B) + Pl''\B) 

+lp^^^\B)]-VMA'^''\B), 
AiB^^c^K*') = -VuV:,[P'''^{B)-\p^^'\B)-'^P^;j'{B) + 2P^*^{B) + Pf'\B) 
-ipX*(5)]. (A-18) 
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For B K*K* decay channels, we have(three part of amplitude are similar): 

A{B'^K*+K*-) = -VuV,l*[Pr^\B)+Pf^\B) + ^Pi^*^\B) 

-^pi§'^\B)] + VudV:,E^'^\B), 

A{B+ ^ K*+K*') = -VtaV,l[P^*^\B) - lpf/*^*(B) + Pf^\B) + 

+V^aV:,A^*^\B), 

A{B'^K*'K*') = -VtaV:,[P^'^\B)-\pq^''^''{B) + Pf^\B) + Pr^\B) 

(A-19) 

Let us first give the factorizablc emission contributions for the {V — A) x {V — A) and 
(y — A) X (y + A) effective four-quark vertexes, they are simply denoted by LL and LR 



TLla^"{M) = i\^FM Fm, Fm, C C C dudxdymlct>M{u) 

4 JVC Jo JO Jo 

{mB(2m6 - mBx)(f)Miix) + i-iMiC^ruBX - m6)[^^(x) - (2;)]}0M2(2/)^ra('"> 
j,FM,M,^^^ = i^^FMFM,FM,j j J dudxdymliXM,(l>M{u)(l>MMKiS^)hTbiu,x), 

T^M,M,(M) = -r£,t^^^(M). (A-20) 

The factorizable emission contributions for the (S — P) x (S + P) effective four-quark vertex are 
found to be 

TsPa'^'iM) = Fm Fm, Fm, C C C dudxdy uib m^^u) 

^ J-^c Jo Jo Jo 

{mB{2mB - mb)(f)Mi {x) + UMi [4mb - {x + l)mB](^^ {x) + HM-^ruB {l-x)(t>lt^{x)} 

<l^M2^y)hTa{u, x), 

Tgp^^^^{M) = Fm Fm-^ Fm2 I / / du dx dy w?b jiM2(t>M{u) 

^ J-^c Jo Jo Jo 

[mBU(f)MAx) + 2(1 - u)iiMi(t^Mii^)WM2iy)^Tb{u,x). (A-21) 
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Similarly, we obtain 

T^^'^^'iM) = \^ Fm Fm, Fm, C C C dudxdymlMn) 

4 ^^c Jq Jq Jo 

{{u- y)mB(l)M,{x) + (1 - x)^Mi[^^Mi(^) + 4>Mi{x)]}(l)M2{y)hTaiu,x,y), 

T^I^'^\M) = -\§^Fm F, Fm, du dx dy mlct>M{u) 

{{u + x + y- 2)mB(t)Mi{x) + (1 - x)fiMi[(t>M^{x) - (lpM^{x)]}(t)M2{y)hTh{u,x,y) 

(A-22) 

for non-factorizable emission contributions with the {V — A) x {V — A) effective four-quark 
vertex, and 



pNMiM, 
-LR 

pNM 
-LRa 



(M) 



1NM1M21 

LRa 
1 Cf 



T'^Rr'\M)+T^m'''\M) 

Fm Fmi Fm2 / / dudxdy m%4>Miu) 
Jo Jo Jo 



rpNMlM2 

^LRb 



(M) 



ANc 

{^Ma miiiiu - X - y + + {u + x-y- lWM^{x)](t^j^^{y) 

-[{u-x-y+ {x) + {u + x-y- l)(f)li^ {x)](j)li^ (y)} 

+{u-y)mB fJ'M2[4^M2iy) - (fM2iyMMAx)]hTa{u,x,y) 

--— - Fm Fmi Fm2 / / dudxdy m%4>M{u) 
4-/Vc Jo Jo Jo 

{/^M2 mi {[{u-x + y)(t)M^ (x) + (n + a; + y - 2)<^^ {x)Wm2 

+[{u -x + y)ct>i,^ {x) + {u + x + y- 2)(/.^^ {x)]4>l,^ (y)} 

+{u + y-l)mB y^M2[0M2(y) + <^M2(y)]<^Mi(a;)|/iT6(«>^>2/) (A-23) 

for non-factorizable emission contributions with the {V — A) x {V + A) effective four-quark 
vertex, and 



pNMiM2 
'-SP 

pNMiM, 
'-SPa 



(M) 
(M) 



■^NMiM2 
SPa 

-TTF' Fmi Fm2 I I I dudxdy m^^MH 
4ivc Jo Jo Jo 



Tg/^iM2 + T^^^^^ (M) 

rl /-l 



{{u + x-y- l)mB(pMi{x) + (1 - x)nMi[(t>M-,{x) - (i)MSx)]}(pM2{y)hTaiu>x,y) 
r^pfi^2(M) = ];^Fm Fm, Fm2 C C C dudxdy ml^M{n) 

4-'Vc Jo Jo Jo 

{{u + y- \)mB4>Mi{x) + (1 - x)hmiWmS^^ + 4>Mi{x)]}4>M2{y)hTh{u, x, ^-24) 



for non-factorizable emission contributions with the {S — P)x (S +P) effective four-quark vertex. 
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We now present the results from annihilation diagram contributions, 

^FMiM2(^) = -I9f Fm Fm, Fm, j I [ dudxdy mlMu) 

4^Vc Jo Jo Jo 

{(1 - y)m%(j)M2{y)4>Mi{x) + 21J.M2 y"Mi[(2 -j/)^5vf2(y) + 2/^M2(y)]^i(^)}^Aa(^>y)> 

^FMiM2(^) = ^^Fm Fm, Fm, C C C dudxdy ml<i>M{u) 
^-'^C Jo JO Jo 

{xm%4>M2{y)<pMi{x) + 2fiM2 /J'Mi[(l + x)(j)lj^{x) - (1 - x)(plj^{x)](plj^{y)}h^^{x,y), 
AFM.M2^^ = ^f^^^(M) + Af^*^^(M), 

^r^r^(M) = Al^,^^^{M) (A-25) 

for the factorizable annihilation contributions with the {V — A)x{y — A) and {V — A)xiy + A) 
effective four-quark vertexes, and 



^FMiM2^^^ = i^M i^Mi Fm2 C C C dudxdy m%ct>M{u) 

^ ^^c Jo Jo Jo 

[(1 - y)iJ-M2WM2iy) + </'M2(y)]<?^Afi(a;) + 2flMi(j)M2{y)(t^Mii^)]h^a{x,y), 

1 

AsPb^^i^) = --^irf Fm, Fm2 / / dudx dy mB(l)Miu){2ij.M2(pMJy)4>Miix) 
^^^c Jo Jo Jo 

+x HMi(l>M2{y)[(l>M^{x) - (pMi{x)]}h^,,{x,y) (A-26) 

for the factorizable annihilation contributions with the {S — P) x {S + P) effective four-quark 
vertex, and 

^^f^^^(M) = A^^^^{M) + K^^'''{M), 

A-LLa^'^'i^'^) = Fm Fmi Fm2 [ [ [ dudxdy m%(pMiu)\[mb + rnBiu-y)]m%4>M2iy)4'Mi{x) 

4 ^Vc Jo Jo Jo ^ 

+MMiAtM2 {[-{u-x-y + l)mB0^^ (x) + {-u -x + y + l)mB(fii, {x)](I>M2 iv)] 
+[(4m6 + {u + x-y- 1)1713) (f^.i^) + {u- x - y + l)mB4>M^{x)](j^j^^{y)^hAa{u,x,y), 

A^u^'iM) = Fm Fm, Fm2 I' I' I' dudxdy ml(t>M{u)[xmlcl>M2i.y)M^) 

+I^Mim2{ -[{u + x + y- + {-u + x-y + l)4>M^{x)](l)M^iy) 

+[{-u + x-y + 1)4>%^ {x) + {u + x + y- 1)4>1^^ {x)WM,{y)}} h^f,{u, x, y) (A-27) 

for the non- factorizable annihilation contributions with the {V — A)x (V — A) effective four-quark 
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vertex, and 



^Ina^'iM) = Fm Fm, Fm, C C C dudxdy mlMu) 

{liM2 [mb + {y- u)mB] [(Am2 (y) - ^Ih ivMMi (x) 

A^^^^^M) = }-^^FMFM,FM,f^ j\udxdy m|c/>M(n){x /xmJ<A^, (x) + (x)]c/>^,^ (y) 
-(1 - n - y)^M2[0M2(y) ~ (t}ii2{y)]<t}hh{x)}hAb{u,x,y) (A-28) 

for the non-factorizable annihilation contributions with the iy—A) x {V —A) and (V—A) x (y+j4) 
effective four-quark vertexes, and 

^SPa^'^i^) = Fm.Fm, ^ ^ dudxdy mBHi{u)[[mb + {x-l)mB]mlct)i,^{y) 

+f^Mim2{'pMiix)[i'^ -x-y+ l)mB(l)\.j^{x) + (-U - X + y + (a^)]'^M2 (y)} 

+ {[4m5 - {-u - X + y + l)mB](i)\.j^{x) - {u - x - y + l)mB(fii^{x)](t)\,^^{y)^hAa{u,x,y), 

<P^*''(M) = Fm Fm, Fm^ C C C dudxdy m|(^M(n)|(-n - y + l)m|0M2(y)<AMi(x) 

4 JVC Jo JO Jo ^ 

+MMi/iAf2 { [(n + X + y - 1)0^^^ (x) - (- U + X - y + l)<f>lj^ {x)](f>M^ (y)} 

+ [(_^ + X - y + l)<Plj^ (x)-iu + x + y- 1)<pIj^ (x)]<,^ x, y) (A-29) 



for the non-factorizable annihilation contributions with the {S — P)x {S+P) effective four-quark 
vertex. 

The functions h\j^ with {A = a, b) from Eqs. to ([OSjl arise from propagators of 

gluon and quark, here Y = F, N denote the factorizable and non-factorizable contributions 
respectively, and X = T,A the emission and annihilation diagrams respectively. They have the 
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following explicit forms: 

1 



[—u[\ — x)m^ — + ie){xm^ — + le) 

[—u[l — x)in^ ~ l^g^ ie)[—um^ — t^-q + ^ej 

"^"^ ' '^■^ {—u{l — x)m'^ — ^i'^ + ie){{l — x){l — u — y)m'^ — fi'^ + ie)' 

'^^^ ' (— u(l — x)m^ — /x^ + ie)((l — — ti)m^ — + ze) ' 

- - + ie)((l - y)m| - ^2 _,_ ' 



X — + ie) 

hAai^,x,y) = (^(1 _ _ ^2 + _ - x)m| - m2 + ie) ' 

hAb(.u,x,y) (-^(--^ _ y^^i _ ^2 _^ -^^(^(^-1^ _ ^ _ y-j^ ^1 _ ^2 (A 30) 

For the P1V2 final states, we can replace wave function in Eqs. ()A-20p to ()A-28p with 

4>P2{x) (t)V2{x), (t>^p^[x) -(t)y^{x), (l)p^{x) (/'y2(x)(2 X - l),/iP2 ^ -771^2 

(A-31) 

in which cPl^ix) = if with = 

For V1P2 final states, we can replace wave function in Eqs. ([X^ to (|X^ with 

4>Pi{x) (t>v-,{x), (pp^x) -(Pv^ix), 4>p^(,x) 0v\(a;)(2 X - l),HPi -mvi , 

m -^J'P2 ■ (A-32) 

For the V1V2 final states, we can replace wave function in Eqs. ()A-20p to ()A-28p with 

4'Pi{x) (pViix), (l)p^ix) -(py^ix), (/)p^(x) 0v\(2;)(2 x - 1),hp^ -my^ , 
(pP^ix) (pViix), 4>%ix) -(Pv^ix), (l^p^ix) (/'y2(x)(2 X - 1),/XP2 my^ . (A-33) 

For the transverse part of two vector meson final states, we list all hadronic matrix elements 
based on six-quark approach(separate to plus/minus parts, which is valuable when considering 
vertex correction) we first get the factorizable emission contributions for the {V — A) x {V — A) 
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{V — A) X {V + ^)effective four-quark vertexes, they are simply denoted by LL and LR 



TlLa+ '(^^) = \^ Pm Fm^ Fm.2 I I I dudx dymBmMimM2(pMiu) 
^-'*C JO Jo Jo 

(2mfe - mBx)(plj^{x)(l)jj,^{y)hfp^{u,x), 
1 Cp 

= oTT" -Pmi i^Mz / / / dudxdymBmM-,mM2(t>M{u) 
^ ^^c Jo Jo Jo 

{(2m6 - mB)mMi(l)MS^^ + mB{2mB - mb)(j)M^{x)}(plj^{y)hTb{u,x), 

Flli^+^{M) = ---^ Fm / / dudxdymBmMimM2<f>M{u) 

^ j-^c Jo Jo Jo 

{mbU(l)'^^{x) - mq(t)M^{x)}(t)Xi^{y)h^^{u,x), 
1 Cp 

^LLb^^^i^) = 7,-f:7- Fm Fmi Fm2 I / / dudx dymBmM^mM2<t^M{u) 
^ ^^c Jo Jo Jo 

{'rn'q4>M^{x) + mB(t>M^{x)}(l)lj^{y)h^b{u, x), 

T[M^M^_^M) = I™ (M), 

Tf^^'^Lm = T^rf^'f'AM). (A-34) 



The factorizable emission contributions for the {S — P) x {S + P) effective four-quark vertex 
are found to be 

sp,+/- y^^-^) ~ ^ spa,+/- v^-^)-^ spb,+/- ^-"^ / ' 
rJ^^;^l(M) = 0, 

T^spl]+%m = 0. (A-35) 
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Similarly, we obtain: 



_ rpNMlM2 



rpNMlM2 

^LLa,-\- 



rpNMlM2 

^LLa,- 



pNMiM2 
'■LLb,+ 



-"■LLb,- 



NM1M2 ^ 

1 /"l /"l 



(M) = 

(M) = 

(M) = 

(M) = 



ICf ^ ^ ^ 

-Fm Fmi 
^^^C Jo Jo Jo 

\%FmFm,Fm2 1 I I 

hTa{u,X,y), 



fff 

Jo Jo Jo 



du dx dymB(pM{u)mMimq(j)lj^ {x)<Pm2 iy)hTai'^^ x, y), 
du dx dym%(j)M{u)mM2{u - y)4>Mi{x)(f>'^^{y) 



du dx dy m%(pM{u)mMi'mM2i—'^ + u + x + y) 



(l>Mi{x)4>M2(y)hTbi'^,x,y) 



dudxdy m\4)M{u){mq(t)ij^{y)[mMAMS^^ 



-2mB(t)li^{x)] - mM24>M2iy)l-'^Mi^Mii^)i-^ + u + x + y) 
+mB{-l + U + y)4>M^{x)]}hTb{u, X, y) 

for non-factorizable emission contributions with the {V — A) x {V — A) effective four-quark 
vertex, and 



(A-36) 



^ LRb,+ y-^^^ I 



_ rpNM\M2 

~ ■^LRa,+/- 



-LRb,+/- 

fl rl rl 



pNMiM2 
LBb- 



(M) 



1 Cp 

T^TT Fmi Fm2 / / dudxdy mBmM2'm'q(t^M{u) 
^J^c Jo Jo Jo 

|0~M2(y)[mMi'/'Mi(3^) - ^B(l>Mi{x)]}h^^{u,x,y) 

■^j;^ Fm Fm^ Fm2 J J J dudx dy mBmMi(t>M{u)^m%{x - l)(l)M^{y) 
+mM2mq4>M^ (y)}<^Mi ix)hTa{u, X, y) 

---^ Fm Fmi Fm2 / / dudxdy mBmM2'm.q(l)M{u) 
2-/Vc Jo Jo Jo 

[(f>~ M2{y)[m,Mi(l)Miix) - mB4>M^ix)]^hTbiu,x,y) 

■^j;^ Fm Fmi Fm2 J J j du dx dy mBmMi(t>M{u)^rn%{x - l)(l)M^{y) 
+mM2^q4'M2 {y)} 4'Mii^)^Tbi^'^'y) (A-37) 



for non-factorizable emission contributions with the {V — A) x {V + A) effective four-quark 
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vertex, and 



rpNMlM2 ( 



_ rpNMlMi 



rpl\lVmVl2 (l\/r\ j_rpNMrM2 ( Tt/f^ 



SPa,+ 



1 Cp 

= -TiTT / / dudx 

^^^c Jo Jo Jo 



(M) 



dy m%(j)Miu)mMimM2 



T™(M) 



1 + u + x- y)(t>M^ {x)(t>M2 (y)^Taiu, X, y) 
1 /.I /.I 



Fm Fm2 / / dudx dy m%(l)M{u){rnq(l)l^^{y)[-mMi(l)Mii^) 
Jo Jo Jo 



rpNMlM2/ Jl/f\ 

^spb,+ y^^-^) 



pNMiM2 
'-SPb- 



(M) 



1 Cf 

+2mB4>M^{x)] + mM2 0M2(y)[-"iMi^!'Mi(^)(-l + u + x-y) 
+mB{u - y)4'M^{x)]}h^^{u,x,y) 
1 Cp 

7:ir^ Fm Fmi Fm2 / / du dx dy mB(l)M{u)m%(l)M{u)'rnMirnq 
2A'c Jo Jo Jo 

(pMiix)<pM2iy)hTbiu,x,y) 

/ dudxdy mB4>M{u)m%4)M{u)mM2 
JO Jo Jo 

(-! + « + {x)^tl, (V)'4i(«, X, y) (A-38) 



for non-factorizable emission contributions with the {S — P)x (S + P) effective four-quark vertex. 
We now present the results from annihilation diagram contributions, 



1 Cp f f f 

---^ Fm Fmi Fm2 / / dudxdy m%mM2(t>Mi: 
^J^ic Jo Jo Jo 

[(f>M2iy)^Mi(^)'^MrX - mg(l)J^^{y)(f)li^{x)]h^^{x,y), 

ICp '•^ 



ICp ^ ^ r. 



„ Fm Fmi Fm2 / / dudxdy m%mMimM2(pM{u) 
2-/Vc Jo Jo Jo 

(l>M2(y^^Mii^)hAaix,y), 

2 / / dudxdy m%niMimM24>Miu) 
Jo Jo Jo 

<^M2(y)^MS^)^Ah{x,y), 

1 Cp 

--— - Fm Fmi Fm2 / / I dudxdy m%mMi(l>M{u) 
^A'c Jo Jo Jo 



(A-39) 

for the factorizable annihilation contributions with the {V — A)x{y — A) and (F — A) x + A) 
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effective four-quark vertexes, and 



5P6,+/- 

1 rl /•! 



Cf 

— irr Fm-, 
J-^c Jo Jo Jo 

Cf 



dudxdy mBTnqmMirnM24'M{u) 



Fm Fmi Fm.2 / / du (h dy mBrnM2(t>Miu)[-<f>M2iy)'^B(t>Mi{x) 

'■^c Jo Jo Jo 



+4>M^{y)'mqmMi(t)M^ (x)]/i^„(x, y), 
-I rl rl 

dudxdy mBmqmMi'mM24>M{u) 



— -TT" Fm Fmi Fm2 

J^c Jo Jo Jo 



llS^)(t>My{x)hAh{X:y) 



Cf f f f 

— Fm Fmi Fm2 / / dudxdy mBmMi(t>M{u)[(t)MAx)m%(t)ii^{y) 
J^^c Jo Jo Jo 

-(t>M2 iy)'rnqmM2(t>Mi ix)]hAbix^ y) (A-40) 



for the factorizable annihilation contributions with the {S — P) x {S + P) effective four-quark 
vertex, and 



AlM.M2^M) 



4TrA(M) + ^-tf/_(M), 

1 Cp 

--— ^ Fm Fmi Fm2 / / dudxdy mBrrn,mM-^mM2(t>M{u)(t)l^ Ax) 
^ J-^c Jo Jo Jo 

(t>M2iy)hAa{u,x,y), 

1 Cf 
~2iVc 



Fm Fmi Fm2 / / dudx dy mBmbmMimM24>M{u)(j)j^ {x) 
Jo Jo Jo 

<I^M2^y)hAa{u,x,y), 
1 Cf 

7,i:r Fm^ Fm2 I / / dudxdy mBmqmM^mM2<t>M{u) 
^J^c Jo Jo Jo 

^tii ix)<pM2{y)hAbiu^ X, y) 

Fm Fmi Fm2 / / / du dx dy mBmqmMimM2<t>M{u) 
Jo Jo Jo 

(A-41) 



1 Cf 

(pM^ix)4>M^iy)hAbiu,x,y) 



for the non-factorizable annihilation contributions with the {V—A)x (V—A) effective four-quark 
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vertex, and 



NMi 
LRa, 



A'lttHM) = 0, 



^LRat'^iM) = \^FmFm,FmJ [ [ dudxdy mlcl^Min) 
^ ^^c Jo Jo Jo 

{mM2 [rrib + mB{-u + y)](t>M.^ iy)4>lh (^) 

+mM^{-mB -mh + mBx)(l)M^{y)(l)lj^{x)}h^^{u,x,y), 



A^M^M^iM) = 0, 



^J^c Jo Jo Jo 

{mMiirriq + mB(-l + n + y)]4>M^{y)4>M^{x) 

-TUMiimq - mBx)4)M^{y)4>l^^{x)}hAb{u,x,y) (A-42) 

for the non-factorizable annihilation contributions with the (V—A) x (V—A) and (V—A) x (y+j4) 
effective four-quark vertexes, and 



A^^lffiM) = <pt!%(M) + <pt*f_(M), 

^SPa + ^i^) = Fm Fm^ Fm2 [ [ [ dudx dy mBinbmMimM^^'Miu) 

^ '-^c Jo Jo Jo 



<pAhi^)(t^M2(y)^Aaiu,X,y), 

"1 /-l rl 



AsPa-'^iM) = Fm Fmi Fm.2 [ [ [ dudx dy mBmbmMirnM2<pM{u) 

^ ^^c Jo Jo Jo 

(I^MSx)<pM2(y)^Aaiu,X,y), 



1 Cp 

A^pI^^'^{M) = --r^ Fm Fm^ Fm2 I I / dudx dy mBmqmMimM2(t>M{u 
^^c Jo Jo Jo 



10 Jo Jo 



<pMii^)(I^M2(y)'^Abiu,x,y) 
1 Cp 

A^p^^^^{M) = Fm Fmi Fm2 / / dudx dy mBmqmMj_mM2(l)M{u) 

^ '■^c Jo Jo Jo 

^liS^)<Pti,{y)h'lb{n,x,y) (A-43) 

for the non-factorizable annihilation contributions with the {S — P)x (S + P) effective four-quark 
vertex. 

The functions /i^^ with {A = a, h) from Eqs. ()A-20p to (|A-28p arise from propagators of gluon 
and quark, here Y = F,N denote the factorizable and non-factorizable contributions respec- 
tively, and X = T,A the emission and annihilation diagrams respectively. Their definitions 
refer to Eqs. (jX30|). 
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